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1.  INTRODUCTION 

This  final  report  covers  the  period  August  30,  1999  through  February  28,  2001.  We  summarize 
our  progress  in  determining  the  delamination  resistance  of  through  thickness  reinforced 
structures  under  dynamic  load  conditions.  Leveraging  collaborations  at  Los  Alamos  National 
Lab  and  the  University  of  Genova,  Italy  and  a  subcontract  with  Dr.  Ares  Rosakis  (ETECH  Inc.) 
for  dynamic  push-in  experiments  will  also  be  described. 

Our  progress  in  this  reporting  period  has  resulted  in  eight  papers  written  or  in  preparation  [1-8]. 
We  are  especially  pleased  with  some  fundamental  results  that  we  have  derived  with  simple 
analytical  models.  The  models  have  yielded  particularly  insightful  results  for  general 
combinations  of  properties  of  the  matrix  and  the  through  thickness  reinforcement,  crack  velocity, 
loading  configurations,  elastic  anisotropy,  and  bridging  parameters. 

Highlights  of  our  accomplishments  are  as  follows. 

(i)  Through-thickness  reinforcement  is  a  promising  solution  to  the  problem  of  delamination 
susceptibility  in  laminated  composites,  but  its  acceptance  by  the  design  community 
awaits  dependable  models  of  its  performance  and  failure  to  treat  mixed  mode  loading 
cases  and  bridging  tows  that  are  canted  relative  to  the  delamination  fracture  plane.  We 
have  constructed  elaborate  micromechanical  models  that  show  how  a  bridging  tow 
should  behave  if  it  is  initially  inclined  to  the  fracture  plane  and  subject  to  mixed  mode 
quasistatic  loading  [1].  From  these  models,  the  effective  bridging  law  for  a  bridged 
delamination  crack  can  be  derived.  Approximations  guided  by  experiments  enable  the 
results  to  be  obtained  in  remarkably  simple  form,  with  closed  form  analytical  expressions 
available  for  certain  cases.  The  approach  outlined  above  is  an  excellent  starting  point  for 
studying  dynamic  delamination  under  mixed  mode  loading  conditions. 

(ii)  The  generalized  model  developed  above  should  be  applicable  to  all  problems  of 
monotonic  loading  and  for  general  stress  state.  The  model  reveals  that,  provided  pullout 
of  the  entire  tow  does  not  occur,  shear  bridging  tractions  are  sustained  more  effectively 
by  fibrous  tows  canted  so  that  they  will  be  loaded  with  the  nap  (tension  along  their  axes) 
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than  by  tows  canted  so  that  they  will  be  loaded  against  the  nap.  Pullout  can  be  avoided 
either  by  using  continuous  stitching  or,  if  short  rods  are  used,  by  ensuring  that  the  aspect 
ratio  of  the  rod  (length/circumference)  is  sufficiently  large.  Orienting  fibrous  tows  to  be 
loaded  in  axial  tension  exploits  the  high  axial  stiffness  of  the  fibres  in  the  tow.  Orienting 
them  so  that  they  are  initially  in  compression  will  cause  initial  fibre  rotation  by  a 
mechanism  analogous  to  kinking  that  depends  primarily  on  the  matrix  strength  and  also 
maximises  the  tendency  of  the  tow  to  deflect  laterally  through  the  laminate.  Both  of 
these  mechanisms  tend  to  create  a  softer  response  under  shear  bridging  tractions. 

A  bridging  law  for  fiber  reinforced  composites  under  dynamic  crack  propagation 
conditions  has  been  derived  [2,3].  Inertial  effects  in  the  mechanism  of  fibre  pullout 
during  dynamic  propagation  of  a  bridged  crack  were  critically  examined  for  the  first  time. 
By  reposing  simple  shear  lag  models  of  pullout  as  problems  of  dynamic  wave 
propagation,  the  effect  of  the  frictional  coupling  between  the  fibres  and  the  matrix  is 
accounted  for  in  a  fairly  straightforward  way.  Analytical  solutions  were  found  to  the 
problem  of  coupled  waves  that  propagate  away  from  the  fracture  plane  of  the  bridged 
crack  as  the  bridging  tractions  increase  with  time.  These  solutions  yield  the  time- 
dependent  relationship  between  the  crack  opening  displacement  and  the  bridging  traction. 

With  the  aid  of  the  dynamic  pullout  law  derived  above,  simple  criteria  were  determined 
for  significant  inertial  effects  in  representative  crack  propagation  problems.  We  find  that 
inertial  effects  in  the  bridging  mechanism  will  often  be  significant  for  a  matrix  crack 
propagating  dynamically  in  the  steady  state  ACK  limit  in  a  brittle  matrix  continuous  fibre 
composite  or  for  delamination  in  laminates  reinforced  by  through-thickness  stitching  or 
rods.  We  also  find  that  significant  inertial  effects  are  favoured  by  low  fibre  volume 
fraction,  low  friction  stress,  low  matrix  bar  wave  speed,  and  low  fibre  modulus;  and  high 
fibre  diameter  and  high  matrix  modulus.  If  the  fibre  modulus  is  high  enough  (relative  to 
the  matrix  modulus),  the  criterion  for  inertial  effects  becomes  independent  of  fibre 
modulus.  The  matrix  density  enters  the  criterion  only  through  the  matrix  bar  wave  speed. 
The  criterion  is  always  independent  of  the  fibre  density.  Inserting  values  for  t5^ical 

stitched  laminates,  one  finds  that  dynamic  effects  will  in  fact  be  significant  for  many 
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delamination  crack  propagation  cases.  Using  a  dynamic  traction  law  rather  than  a  law 
for  static  loading  is  therefore  essential  in  general. 

(v)  For  pullout  or  bridging  stresses  that  linearly  increase  in  time,  the  instantaneous  crack 
displacement  is  less  in  the  presence  of  inertial  effects  than  it  would  be  under  static 
loading  to  the  same  bridging  stress.  However,  solutions  for  pullout  from  a  rigid  matrix 
suggest  that,  if  the  bridging  stress  rises  rapidly  and  is  then  held  at  a  constant  value,  the 
crack  displacement  when  all  particle  motion  finally  stops  will  be  greater  than  it  would 
have  been  under  static  loading  to  that  stress  level.  Thus  regimes  of  both  hardening  and 
softening  of  the  bridging  traction  law  due  to  inertial  effects  can  be  expected  in  bridged 
crack  problems. 

(vi)  We  also  formulated  fast  numerical  solutions  for  general  loading  cases  and  hence  have  a 
scheme  to  compute  the  bridging  law  for  fairly  general  conditions  [4].  This  will  be 
required  to  solve  large-scale  bridging,  dynamic  crack  problems  to  self-consistency 
(bridging  tractions  unknown  a  priori).  The  history  dependent  bridging  law  thus 
computed  will  form  the  basic  ingredient  in  the  computation  of  the  structural  response  of 
the  composite  under  dynamic  loading  conditions. 

(vii)  The  dynamic  delamination  cracking  behavior  and  the  energetics  of  crack  grovrth  in 
through  thickness  double  cantilever  beam  (DCB)  specimens  has  been  analyzed  [5].  The 
double  cantilever  beam  (DCB)  specimen  loaded  dynamically  by  a  flying  wedge  offers  a 
relatively  simple  experimental  approach  to  analyzing  the  mode  I  dynamic  delamination 
problem.  The  role  of  bridging  by  stitches  or  rods  on  dynamic  crack  growth  was  computed 
by  solving  the  bridged  crack  problem  within  the  framework  of  beam  theory.  Analytical 
results  were  obtained  for  steady  state  crack  propagation  conditions. 

(viii)  The  magnitude  of  the  bridging  traction,  which  is  controlled  by  the  properties  of  the 
through  thickness  reinforcement,  plays  a  significant  role  in  the  crack  growth  behavior. 
For  steady  state  crack  growth  conditions,  different  regimes  of  the  solution  behavior  have 
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been  identified  which  would  correspond  to  different  crack  propagation  characteristics. 
Critical  velocities  that  demarcate  the  various  solution  regimes  have  been  identified. 

(ix)  Steady  state  crack  growth  is  attainable  for  this  loading  configuration  provided  certain 
conditions  are  satisfied.  Regions  of  stable  crack  growth  have  been  identified  in  terms  of 
the  material  properties  of  the  through  thickness  reinforcement,  the  size  of  the  DCB 
specimen  and  the  velocity  of  the  wedge.  This  provides  guidelines  for  design  of 
experiments  to  probe  the  efficacy  of  bridging  on  improving  the  dynamic  fracture 
toughness  of  through  thickness  reinforced  structures  [6]. 

(x)  The  dynamic  crack  energy  release  rate  as  a  function  of  the  local  crack  tip  velocity  has 
been  computed  for  a  material  possessing  orthotropic  symmetry  under  general  mixed 

mode  loading  conditions.  The  dynamic  energy  release  rate  varies  as  0[(v,,  —  v)”*]as 
V  V,. ,  and  where  v  is  the  crack  speed  and  Vr  is  the  Rayleigh  wave  speed.  We  find  that 
the  variation  in  the  dynamic  crack  energy  relase  rate  for  orthotropic  materials  can  be 
rationalized  in  terms  of  the  orthotropy  parameters  X  and  p,  just  as  in  the  static  case.  This 
computation  of  the  dynamic  energy  release  rate  is  necessary  to  calculate  the  energetics  of 
crack  growth  for  standard  engineering  specimens  [7]. 

(xi)  An  approximate  formula  for  the  static  weight  function  that  describes  the  crack  tip  stress 
intensity  factor  due  to  a  pair  of  point  loads  on  the  fracture  surface  of  an  orthotropic 
delamination  specimen  has  been  postulated  and  validated  numerically  [8].  This  weight 
function  enables  accurate  solution  of  mode  I  delamination  problems  in  the  presence  of 
large  scale  bridging  over  all  crack  length  regimes,  since  it  avoids  the  errors  of  beam 
theory  in  describing  the  singular  crack  tip  fields.  However,  comparison  of  results  for 
large  scale  bridging  problems  with  a  Dugdale  (uniform  traction)  bridging  law  shows  that 
a  beam  theory  model  with  appropriate  crack  tip  corrections  (available  in  the  literature) 
may  be  acceptably  accurate.  The  weight  functions  will  be  used  to  formulate  integral 
equation  solutions,  which  will  be  to  ensure  that  accurate  results  are  always  available 
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least  as  a  reference  in  future  work.  Extension  to  mode  II  loading  and  dynamic  weight 
functions  is  under  way. 

(xii)  A  subcontract  was  issued  to  Etech,  Inc.  for  experiments  to  gain  understanding  of  the 
mechanics  of  fibre  push-in  (analogous  to  pullout)  under  dynamic  loading.  The  first 
experiments  have  already  been  completed  (Owen  and  Rosakis)  and  the  key  strain  field 
and  velocity  data  successfully  recorded.  Further  reports  will  follow  after  proper  analysis 
of  the  data. 

Considerable  headway  has  been  made  in  all  of  the  tasks  listed  above,  even  though  the  program 
has  been  in  effect  for  only  a  short  time  period  (~  16  months).  The  work  in  each  case  is  mainly 
finished  and  we  are  well  advanced  in  writing  papers.  Unfinished  mansucripts  will  be  completed 
in  the  next  few  months. 
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2.0  FUTURE  WORK 

The  following  issues  remain  as  the  natural  continuation  of  our  work. 

1.  The  dynamic  traction  law  model  must  be  generalised  to  mixed  mode  loading.  We  have  a 
model  for  mixed  mode  loading  in  the  static  case  [1].  We  now  know  from  this  year’s 
work  that  dynamic  effects  are  likely  to  be  important  in  most  cases  and  we  therefore 
should  consider  how  to  adapt  our  static  model  to  deal  with  inertial  effects  during  mode  II 
deflections  of  a  stitch  in  a  laminate. 

2.  The  fracture  or  crack  propagation  problem  must  be  solved  for  mixed  mode  loading.  We 
know  from  the  static  case  that  mode  ratio  effects  have  a  major  influence  on  the 
qualitative  character  of  crack  propagation  in  the  presence  of  large  scale  bridging.  This 
fundamental  question  must  therefore  also  be  resolved  for  dynamic  fracture. 

j.  We  are  now  in  position  to  map  out  (a)  a  method  of  predicting  the  extent  of  delamination 
damage  during  composite  impact  and  (b)  the  energy  absorption  capacity  of  delaminating 
composites.  In  both  cases,  the  outcome  will  be  very  strongly  affected  by  the  presence  of 
through-thickness  reinforcement  (large  scale  bridging  zones). 

4.  Experimental  data  must  be  obtained  and  analyzed  for  dynamic  crack  and  pullout 
problems  to  validate  our  theoretical  results.  Although  there  has  been  considerable 
experimental  characterization  of  fibre  -  matrix  debonding  under  quasi-static  loading 
conditions,  comparatively  little  is  known  regarding  the  process  at  relatively  high, 
dynamic  loading  rates.  Prof  Ares  Rosakis  of  Caltech,  with  a  sub-contract  under  this 
program,  is  conducting  a  detailed  experimental  study  of  the  dynamic  fiber  pull-out  / 
push-in  process  and  will  provide  data  with  which  the  validity  of  the  pull-out  model  will 
be  assessed.  Through  the  use  of  a  carefully  designed  2D  model  material  system  and 
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high-speed  diagnostics,  a  complete  quantitative  description  of  process  can  be  obtained 
which  makes  direct  comparison  with  the  model  possible. 

5.  Upon  passing  of  the  crack  tip,  a  debond  crack  propagates  along  the  length  of  the  fibre 
away  from  the  fracture  plane.  Propagation  of  the  debond  crack  is  governed  by  the 
fracture  energy  associated  with  the  separation  of  the  matrix  and  the  fibre  at  the  debond 
crack  tip  and  the  work  done  against  friction  in  displacing  the  debonded  fibre  along  its 
axis.  In  many  composites,  the  debond  energy  is  small  and  pullout  is  dominated  by 
friction  over  a  wide  range  of  pullout  displacements.  However,  to  deal  with  cases  where 
this  assumption  is  not  true,  we  intend  to  develop  a  theoretical  framework  that 
incorporates  the  debond  energy  in  the  pullout  model.  Although,  simple  extensions  of 
shear  lag  models  are  available  to  deal  with  debond  energy  contribution  for  the  quasistatic 
loading  case,  the  corresponding  formalism  for  the  dynamic  case  appears  much  more 
complicated  and  is  yet  to  be  developed. 

6.  Weight  function  or  other  computational  formulations  must  be  used  to  validate  and  extend 
anal5^ical  and  illustrative  results  based  on  beam  theory. 

7.  Based  on  our  increasing  understanding  of  the  dynamic  delamination  process  with  large 
scale  bridging,  models  of  energy  absorption  in  a  laminate  structure  and  delamination 
resistance  must  be  formulated.  Such  models  are  intended  to  enable  systematic 
optimisation  of  energy  absorption  and  delamination  resistance  by  providing  reasonably 
simple  but  physically  correct  design  rules. 

These  are  the  issues  identified  in  our  original  proposal.  We  are  very  much  on  track  to  address 

them  in  their  turn. 
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3.0  SUMMARY  OF  PROGRESS 


3.1  Mechanics  of  Crack  Bridging  by  Fiber  Tows  under  Mixed  Mode 
Loading 

Through-thickness  reinforcement  is  a  promising  solution  to  the  problem  of  delamination 
susceptibility  in  laminated  composites,  but  its  acceptance  by  the  design  community  awaits 
dependable  models  of  its  performance  and  failure.  Models  have  been  slow  in  coming,  mainly 
because  the  mechanisms  involved  in  the  failure  of  through-thickness  reinforcing  elements  are  at 
first  sight  very  complicated. 


However,  from  a  micromechanical  standpoint,  key  mechanisms  of  damage  recently  revealed  by 
detailed  experiments  show  encouraging  universality  among  bridging  tows  of  various  kinds, 
including  stitches  and  short  fibrous  rods  (so-called  “z-fibers”).  The  key  phenomena  are  as 
follows:  The  bridging  tow  deforms  plastically  in  axial  shear  at  the  smallest  observable  crack 
sliding  (mode  II)  displacements.  Plasticity  is  mediated  at  first  by  crazing  and  microcracking  in 
the  resin  within  the  tow.  Plasticity  and  damage  (microcracking,  splitting,  and  spalling)  occur 
soon  after  in  the  adjacent  laminate.  Plasticity  allows  fiber  rotation  near  the  delamination  crack 
plane,  so  that  the  axial  stress  in  the  fibers  creates  significant  shear  components  on  the  fracture 
plane  (creating  shear  bridging  tractions).  By  this  mechanism,  the  mode  II  bridging  tractions  can 
reach  ~  1  GPa  across  the  section  of  the  tow,  which  is  an  order  of  magnitude  higher  than  the  tow 
could  support  in  axial  shear  prior  to  fiber  rotation.  At  small  crack  opening  or  sliding 
displacements,  the  tow  also  detaches  via  circumscribing  matrix  cracks  (not  necessarily  fiber/resin 
debond  cracks)  from  the  surrounding  composite.  The  debonded  tow  slides  through  the  laminate 
and  the  debonded  zone  extends  fairly  quickly  to  the  outer  surface  of  the  laminate.  Sliding 
displacement  can  accommodate  the  increase  in  path  length  that  occurs  when  the  locus  of  the 
bridging  tow  is  altered  by  fiber  rotation.  Fiber  rotation  in  the  bridging  tow  also  implies  lateral 
motion  of  the  tow  into  the  surrounding  laminate.  This  ploughing  action  may  be  modeled 

phenomenologically  as  the  problem  of  a  punch  being  driven  into  a  plastic  medium,  since  the 
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laminate  deforms  by  matrix  mediated  plasticity.  The  driving  force  for  the  punch  deflection  may 
be  regarded  as  the  combination  of  the  line  tension  sustained  by  the  fibers  in  the  bridging  tow  in 
concert  with  their  curvature,  which  increases  as  the  tow  deflects. 

In  this  work,  all  of  the  phenomena  described  above  has  been  modeled  with  remarkable 
simplicity,  and  without  betraying  the  observed  physics  of  the  bridging  process  [1].  A  model  of 
the  response  of  the  bridging  tow  has  been  formulated  and  solved  analytically,  giving  useful 
insight  into  the  governing  material  and  geometrical  factors.  The  formulation  is  generalised  to 
allow  the  bridging  tow  to  be  canted  relative  to  the  fracture  plane  and  to  deal  with  mixed  mode 
loading  conditions.  The  essential  material  property  that  came  out  of  solving  the  mechanics  of 
the  bridged  crack  problem  is  the  traction  law,  p(u),  relating  the  vector  of  relative  displacements 
of  the  crack  surfaces,  u,  to  the  vector  of  tractions,  p,  supplied  to  the  fracture  surfaces  by  the 
bridging  through-thickness  reinforcement.  Generalizations  to  any  combination  of  mode  I  and  II 
displacements  and  tractions  was  considered  to  make  it  relevant  to  problems  of  interest.  The 
mixed  mode  traction  law  was  derived  in  analytical  form,  requiring  at  most  the  solution  of  two 
transcendental  equations  in  a  single  scalar  variable.  The  results  contain  only  a  few  material 
parameters,  including  a  shear  flow  stress  for  the  bridging  tow,  a  friction  stress  for  tow  sliding, 
and  a  punch  resistance  stress,  all  of  which  can  be  independently  estimated.  Reasonable  values  for 
the  parameters  lead  to  laws  that  agree  very  well  with  experimental  measurements.  Thus  a  good 
description  is  readily  available  for  the  physics  of  crack  bridging  in  the  static  delamination  crack 
problem. 

The  generalized  model  developed  above  should  be  applicable  to  all  problems  of  monotonic 
loading  and  for  general  stress  state.  The  model  reveals  that,  provided  pullout  of  the  entire  tow 
does  not  occur,  shear  bridging  tractions  are  sustained  more  effectively  by  fibrous  tows  canted 
such  that  the  tows  are  loaded  with  the  nap  (tension  along  their  axes).  Pullout  can  be  avoided 
either  by  using  continuous  stitching  or,  if  short  rods  are  used,  by  ensuring  that  the  aspect  ratio  of 
the  rod  (length/circumference)  is  sufficiently  large.  Orienting  fibrous  tows  to  be  loaded  in  axial 
tension  exploits  the  high  axial  stiffness  of  the  fibres  in  the  tow.  Orienting  them  so  that  they  are 
initially  in  compression  will  cause  initial  fibre  rotation  by  a  mechanism  analogous  to  kinking  that 
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depends  primarily  on  the  matrix  strength  and  also  maximises  the  tendency  of  the  tow  to  deflect 
laterally  through  the  laminate.  Both  of  these  mechanisms  tend  to  create  a  softer  response  under 
shear  bridging  tractions. 

To  conclude,  we  have  constructed  an  elaborate  micromechanical  model  that  shows  how  a 
bridging  tow  should  behave  if  it  is  initially  inclined  to  the  fracture  plane  and  subject  to  mixed 
mode  loading.  From  these  models,  the  effective  bridging  law  for  a  bridged  crack  model  of  the 
delamination  can  be  derived.  Once  again,  approximations  guided  by  experiments  enable  the 
results  to  be  obtained  in  remarkably  simple  form,  with  closed  form  anal)^ical  expressions 
available  for  certain  cases.  The  approach  outlined  above  is  an  excellent  starting  point  for 
studying  dynamic  delamination  under  mixed  mode  loading  conditions. 


3.2  Inertial  effects  in  the  pullout  mechanism  of  a  bridged  crack 

In  this  part  of  our  research  [2,3],  it  is  our  aim  to  develop  an  equivalent  bridging  model  for 
individual  through-thickness  reinforcement  tows  under  dynamic  loading.  Dynamic  loading  rates 
can  significantly  affect  the  mechanisms  described  above,  probably  leading  to  a  hardening  of  the 
traction  law  and  more  brittle  behavior  of  the  bridging  tow.  In  order  to  get  insights  into  how 
dynamic  effects  modify  the  mechanics  of  deformation,  we  begin  by  modifying  the  mechanics  of 
pullout  under  quasistatic  loading  to  include  the  inertial  effects. 

The  mechanics  of  pullout  and  the  resulting  traction  law  have  been  much  studied  and  are  well 
understood  for  static  loading.  Simple  analytical  forms  are  available  for  p(ti)  when  the  frictional 
coupling  of  the  reinforcement  to  the  matrix  is  uniform  and  slip  extends  over  distances  that  are 
large  compared  to  the  reinforcement  diameter  (Marshall,  Cox,  and  Evans,  1985;  McCartney, 
1987).  In  this  limit,  which  is  a  common  case  in  ceramic  composites  and  textile  polymeric 
composites,  the  shear  lag  model  of  load  transfer  between  the  reinforcement  and  the  matrix  is 
accurate.  Simple  extensions  of  models  of  this  class  are  also  available  to  deal  with  small  but 
nonzero  levels  of  the  work  required  to  debond  the  reinforcement  from  the  matrix  prior  to  slip 
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[11],  Given  the  relationship,  p{u),  the  characteristics  of  crack  propagation  can  be  computed  by 
solving  a  bridged  crack  problem  [19-22]. 

In  this  section,  we  extend  existing  models  of  the  mechanics  of  pullout  to  high  loading  rates.  An 
approach  to  evaluating  a  traction  law  that  takes  account  of  the  inertia  of  the  reinforcement  and 
the  matrix  is  formulated  as  a  direct  extension  of  the  elementary  static  loading  model  [9,10].  The 
chosen  base  model  for  static  loading  has  proven  consistent  with  experiments  in  many  material 
systems  and  has  been  the  foundation  of  major  advances  in  understanding  damage  in  composites. 
The  spirit  of  the  present  work  is  to  seek  equivalent  insight  into  dynamic  damage  by  incorporating 
the  influence  of  inertia  into  the  simplest  credible  model.  Thus  not  all  aspects  of  the 
micromechanics  of  pullout  that  could  be  important  in  some  cases  will  be  addressed.  Instead, 
attention  will  be  focused  on  identifying  a  characteristic  time  for  the  frictional  pullout  problem 
that  will  allow  rapid  assessment  of  when  inertial  effects  will  be  important. 

Idealization  of  the  Bridged  Crack  problem 

The  crack  propagation  and  pullout  problems  are  depicted  schematically  in  Fig.  1.  A  matrix 
crack  propagates  on  the  plane  z  =  0  and  is  bridged  by  intact  fibres  in  its  wake.  (For  simplicity  of 
expression,  the  term  fibre  from  here  on  will  be  used  to  refer  to  bridging  entities  of  any  kind, 
including  stitches,  rods,  and  bridging  grains.)  Upon  the  passing  of  the  crack  tip,  a  debond  crack 
propagates  along  the  length  of  each  fibre  away  from  the  fracture  plane.  Propagation  of  the 
debond  crack  is  governed  by  the  fracture  energy  associated  with  the  separation  of  the  matrix  and 
the  fibre  at  the  debond  crack  tip  and  the  work  done  against  friction  in  displacing  the  debonded 
fibre  along  its  axis.  In  many  composites,  the  debond  energy  is  small  and  pullout  is  dominated  by 
friction  over  much  of  the  range  of  pullout  displacements.  Therefore,  in  this  first  model  the 
debond  energy  will  be  ignored. 

For  static  loading,  the  crack  propagation  problem  can  be  idealised  by  replacing  the  process  zone 
by  elastic  composite  material  down  to  the  fracture  plane  and  representing  the  phenomena  within 
the  process  zone  by  bridging  tractions  applied  continuously  on  the  fracture  surfaces  (Fig.  1).  The 
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bridging  tractions,/?,  are  related  to  the  axial  stress,  T,  in  the  fibres  at  the  fracture  plane  by p  =fr, 
where  f  is  the  area  fraction  of  the  fibres  on  the  plane  z  —  0.  For  aligned  continuous 
reinforcement, /is  also  the  volume  fraction  of  the  fibres.  The  total  crack  opening  displacement, 
2u,  in  the  idealisation  should  be  defined  as  the  difference  in  the  actual  displacement  evaluated 
across  the  process  zone  and  the  displacement  that  would  be  expected  if  the  material  in  the 
process  zone  was  elastic. 
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Figure  1.  (a)  Schematic  of  a  crack  bridged  by  fibres,  showing  process  zone  where  relative 

displacement  exists  between  fibres  and  matrix,  (b)  Idealization  of  the  bridged  crack  problem 
with  process  zone  replaced  by  surface  tractions  acting  on  the  fracture  surfaces. 

The  traction  law,  p(u),  is  derived  by  considering  the  micromechanics  of  the  phenomena 
occurring  within  the  process  zone,  which  is  to  say  the  micromechanics  of  frictional  sliding.  The 
micromechanical  problem  can  be  represented  by  a  small  volume  of  material,  e.g.,  the  materia! 
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bounded  by  one  of  the  dotted  rectangles  in  Fig.  la.  The  traction  boundary  conditions  for  the 
representative  volume  are  as  follows.  At  2  =  0,  the  matrix  is  traction-free,  while  the  fibres 
sustain  the  axial  traction  T.  At  2  =  4,  the  strain  in  the  fibres  and  the  matrix  must  equal  the 
average  strain  in  the  composite  adjacent  to  the  process  zone  (2  >  4).  Shear  tractions  may  arise 
along  the  vertical  boundaries  of  the  representative  volume  (parallel  to  2),  but  these  are  neglected. 

In  the  problem  of  a  crack  propagating  under  static  loading,  the  traction  at  z  =  4  for  the 
representative  volume  associated  with  a  fixed  volume  of  material  is  assumed  to  rise  from  zero 
when  the  material  is  immediately  behind  the  crack  tip  to  increasing  values  as  the  material  passes 
further  back  into  the  wake  [9].  The  tractions  are  not  truly  zero  right  at  the  crack  tip,  since  the 
fibres  are  not  stress-free,  but  assuming  they  are  leads  to  reasonable  results  for  the  bridged  crack 
problem  as  long  as  the  fibres  remain  intact  over  sufficiently  long  distances  into  the  crack  wake. 
Then  shielding  of  the  crack  tip  is  dominated  by  the  bridging  tractions  acting  in  the  further  crack 
wake,  where  the  boundary  conditions  are  correct.  In  dynamic  loading,  the  boundary  conditions 
at  the  boundary  of  the  process  zone  involve  displacement  and  displacement  rates  as  well  as  stress 
or  strain  conditions. 

The  Micromechanics  of  Dynamic  Pullout 

When  friction  is  the  only  active  force  of  resistance,  the  static  pullout  problem  reduces  to  an 
idealization  in  which  fibers  are  pulled  out  of  a  half-space  to  which  they  are  not  initially  bonded. 
The  dynamic  problem  is  one  of  wave  propagation  along  a  fibre  in  the  presence  of  frictional 
retardation  and  with  no  debond  crack  tip  or  crack  tip  field  to  be  considered. 

The  archetypal  problem  is  illustrated  in  the  figure  2  below.  A  representative  volume  consists  of 
cylindrical  fibres  of  radius  R  and  volume  fraction  f  embedded  in  a  matrix  (2  >  0).  The  fibre  and 
the  matrix  have  axial  Young’s  moduli  E{  and  Em  and  densities  yOf  and  pm  respectively.  The  axial 
displacement,  strain,  and  stress  of  the  fibre  and  the  matrix  are  denoted  U{  and  Um,  Sf  and  Sm,  and 
Of  and  (Jm,  respectively.  The  axial  displacements  will  be  assumed  to  be  the  only  nonzero 
displacement  components  induced  by  loading  and  to  be  uniform  across  any  section  of  the  fibre  or 


13 


Esciciiiieii 
Science  Center 

SC71164.RPRTTA 

the  matrix.  Thus  the  displacement,  strain,  and  stress  in  the  fibre  and  the  matrix  are  functions  of  z 
and  /  only.  They  are  consistent  with  assuming  that  the  friction  forces  are  constant  (unaffected  by 
fibre  contraction  or  dilation  due  to  axial  stresses). 


7' 

matrix  crack 
plane 


Rscosae  coxoeosoo 


Figure  2:  Schematic  of  the  dynamic  pullout  problem  in  a  composite  near  the  fracture  plane. 


There  is  no  initial  bond  between  the  fibres  and  the  matrix.  The  fibres  are  coupled  to  the  matrix 
by  friction  tractions,  n,  which  are  assigned  the  following  properties.  If  relative  motion  exists 
between  the  fibres  and  the  matrix,  then 

rf=r  (3.2.1a) 

rf=-r(Wf>M„)  (3.2.1b) 

where  a  dot  indicates  time  differentiation,  r  is  a  positive  constant  and  rf  >  0  indicates  friction 
tractions  acting  on  the  fibres  in  the  positive  z  direction.  When  the  fibres  and  the  matrix  are  not  in 
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relative  motion,  the  friction  tractions  may  support  stress  gradients  in  the  fibres  and  the  matrix, 
provided  that  the  required  magnitude  of  rf  does  not  exceed  r.  Thus,  by  equilibrium 
considerations. 


5(Tf  diif 
Pf  — 


dz 


dt 
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dz 


dt 


2/ 


\-fR 


(3.2.1c) 


Here  the  possibility  is  included  that  the  matrix  and  the  fibres  have  the  same  non-zero  velocity 
and  may  also  be  accelerating  together,  although  such  general  solutions  will  not  be  exhibited  in 
this  paper.  With  such  a  friction  law,  the  dynamic  wave  equations  describing  those  parts  of  the 
fibres  and  the  matrix  that  are  in  relative  motion  may  be  written  approximately  as 


d^Ui  _  _  26r  1  d^Uf 

dz^  RE^  Cf  dt^ 


(3.2.2a) 


and 


_  2/  ^  1  d\i^ 

dz^  \-fRE^  cl  dt^ 


(3.2.2b) 


where  ^  -  1  if  u^>  i/^and  9=  ~\  otherwise;  and  C{  and  Cm  are  the  bar  wave  velocities  in  the 
fibres  and  the  matrix,  given  by 


Cf 
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(3.2.2c) 
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The  approximation  of  using  the  bar  wave  velocities  rather  than  the  longitudinal  wave  velocities 
in  the  z  direction  in  the  wave  equations  is  consistent  with  the  simplified  treatment  of  stresses  and 
strains  and  the  assumed  uniformity  of  the  friction  stress. 

Boundary  conditions  in  the  dynamic  case  are  as  follows.  At  the  fracture  plane, 

?/f(0, 0  =  0  &  am(0, 0  =  0  (z  =  0).  (3.2.3) 

At  the  boundary  of  the  process  zone  (limit  of  relative  fibre/matrix  motion), 

«f  =  Wm  =  tic  (z  =  Is)  (3.2.4) 

where  Hc  is  the  displacement  of  the  adjacent  intact  composite;  and  conditions  also  exist  on  stress 
or  strain  and  particle  velocities.  These  further  conditions  depend  on  the  nature  of  the  loading 
history,  which  can  be  expressed  as  the  function  s(t),  where  sis  the  strain  in  the  z  direction  in  the 

intact  composite  adjacent  to  the  process  zone  boundary.  The  bridging  traction,  p,  is  related  to  s 
by 


P  =  ^  (3.2.5a) 

where  the  composite  modulus,  E,  is  given  by 

E=fEi +(\-f)Era  (3.2.5b) 

In  the  depiction  of  Fig.  2,  the  process  zone  boundary  will  propagate  away  from  the  fracture  plane 
as  s(t)  rises.  In  the  case  to  be  considered  in  this  paper,  dj)  will  be  assumed  to  rise  continuously 
and  monotonically  from  zero.  Placing  the  origin  of  time,  /  =  0,  at  the  onset  of  nonzero  e/if),  the 
location  of  the  zone  boundary  at  time  t  may  then  be  written  as: 

h=ri(f)Cta  (3.2.62) 
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where  the  function  T]{f)  depends  on  material  and  geometrical  parameters  and  the  form  of  sit)  and 
will  be  shown  to  be  bounded  by 

0  <  77(0  ^  1  •  (3.2.6b) 

For  the  loading  conditions  considered,  the  additional  boundary  conditions  at  z  =  4  are: 

7-,f=,?m=f/c  (2  =  4)  (3.2.7a) 

and 

et^e^  =  sit)  (z  =  /s).  (3.2.7b) 

If  the  load  history,  sit),  possesses  discontinuities,  e.g.,  a  step  load,  then  discontinuities  in  stress 
and  velocity  will  also  propagate  at  the  boundary  of  the  process  zone  (e.g.,  Achenbach  [12] ). 

Composite  Stress  Risine  Linearly  in  Time 

A  case  of  representative  interest  for  dynamic  bridged  crack  problems  and  for  which  analytical 
results  can  be  found  is  that  of  a  load  or  bridging  traction  that  increases  linearly  in  time.  A 
linearly  increasing  load  might  give  insight,  for  example,  into  bridging  effects  in  a  specimen  in 
which  substantial  bending  arises,  such  as  a  standard  double  cantilever  beam  delamination 
specimen.  In  such  specimens,  the  crack  profile  is  often  approximately  linear  and  the  rate  of 
increase  of  the  bridging  tractions  at  any  point  might  also  be  approximately  linear  if  the  crack 
propagates  at  approximately  constant  speed.  The  bridging  traction  at  a  particular  material  point 
might  be  expected  to  rise  from  zero  as  the  delamination  crack  first  passes  until  a  peak  value  is 
reached,  perhaps  corresponding  to  bridging  fibre  rupture.  Analytical  results  can  be  found  for 
linearly  increasing  loads. 
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e(t)=kt 


(3.2.8) 


where  k  is  constant  and  all  displacements  and  boundary  tractions  are  zero  for  r  <  0.  In  this  case, 
9=  1  and  the  wave  equations  have  the  solutions 


Uf=-—ccz^+[l  +  2Tja]£z 


(3.2.9a) 


(3.2.9b) 


where  t]  is  independent  of  time  and  satisfies 


aif  +  Yzif  +  Par}  -1/2  =  0 


with  the  dimensionless  parameters  a  and  yff  given  by 

a  =  ^  &  P  =  ^El. 

l-fE^ 

Analysis  shows  that  Eq.  (3.2.9c)  has  only  one  real  root,  which  always  lies  in  (0,1). 


(3.2.9c) 


(3.2.9d) 


The  particle  velocities  and  the  strains  for  0  <  z  <  4  are  given  by 


Wf  =  [l  +  2T]a^z 


(3.2.10a) 


"m  =[l  +  2r7aK 


(3.2.10b) 


=  1  +  2r]a\  ^ 


(3.2.10c) 
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s„=ys  (3.2.  lOd) 

The  matrix  particle  velocity  is  uniform  and  increases  in  proportion  to  U.  The  fibre  velocity  rises 
linearly  through  the  process  zone.  A  strain  concentration,  £f/s,  propagates  in  the  fibre  behind  the 
process  zone  front,  while  the  matrix  strain  in  the  process  zone  is  always  less  than  e.  The  strain 
distributions  and  velocities  beyond  the  process  zone,  z  >  4,  need  not  be  specified,  provided  the 
composite  strain  £  obeys  the  condition  Eq.  (3.2.7)  on  the  process  zone  boundary. 

The  displacement,  w,  to  be  used  in  defining  the  traction  law,  is  given  by  the  common  fibre 
and  matrix  displacement,  ui  =  iif(z  =  4)  =  Um{z  =  4)  at  the  boundary  of  the  process  zone  minus 
the  displacement  expected  if  the  process  zone  material  were  elastic: 

u  =  7/1  -  els 


(3.2.11) 


Equation  (3.2.11)  constitutes  the  traction  law  for  the  case  of  linearly  rising  loads.  The  limit  of 
very  fast  loading  corresponds  to  4  — >  oo,  whereupon  a  — >  0  and  77  ->  1,  since  the  first  and  third 
terms  in  Eq.  (3.2.9c)  become  negligible.  The  disturbance  then  propagates  at  the  bar  wave  speed 
in  the  matrix.  Static  loading  is  represented  by  the  limit  4  ->  0  or  a  -»  00,  for  which  the  first  two 
terms  of  Eq.  (3.2.9c)  become  small  and  one  has  the  asymptotic  solution 


7 


(3.2.12) 


For  this  limit,  substituting  Eq.  (3.2.12)  into  Eqs.  (3.2.6a)  and  (3.2.1 1)  yields 
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/  2t 


(3.2.13a) 


/  J  4E,r 


(3.2.13b) 


which  coincide  with  the  results  obtained  by  McCartney  [10]. 

Equations  (3.2.10)  and  (3.2.13b)  show  that  for  loads  that  increase  linearly  with  time,  the  form  of 
the  traction  law  is  identical  in  the  static  and  dynamic  cases  and 


■  =  2afirj 


(3.2.14) 


Thus  the  strength  of  the  effects  arising  from  inertia  is  measured  by  the  degree  to  which  the  limit 
of  Eq.  (12)  is  not  approached.  This  can  be  conveniently  summarised  by  comparing  the  product 
lafit]  to  unity. 

Large  p.  Asymptotic  analysis  of  Eq.  (3.2.9c)  shows  that  in  the  limit  oo,  laPrj  -> 


'^{ccPY  +{apy'  -{apY  (second,  third,  and  fourth  terms  of  Eq.  (3.2.9c)  dominant),  which  is 

a  function  of  ap  and  not  of  ^  separately;  and  Fig.  3a  shows  that  this  limit  is  approached  quite 
closely  for  1,  which  is  expected,  for  example,  for  composites  containing  relatively  stiff  fibres 
(Ef  >  E^m).  Fig.  3  a  also  suggests  that,  as  an  engineering  estimate,  inertial  effects  are  large  when 


ap  =  ~ — <  2 


{P>  1) 


Remembering  that  k  is  the  time  constant  of  the  loading,  this  condition  can  be  rewritten 


f  T  C 


(3.2.15) 


(3.2.16) 
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(a) 


ascosajcooowwj 


(b) 


R$O»»COX06C'7lli 


Figure  3:  The  product  2aPri,  which  indicates  the  relative  importance  of  inertial  effects  in  the 
bridging  phenomenon,  (a)  Numerical  results  for  representative  values  of  /?.  (b)  Numerical 
results  for  laprf  compared  with  asymptotic  limit  for  large  Era  (dashed  curves).  The  value  of 
{aP)\  is  marked  for  P=  0.01  and  is  to  be  compared  with  the  value  of  aP  at  which  laPt] «  0.9, 
which  is  taken  as  a  representative  cutoff  for  significant  inertial  effects. 
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Small  p.  The  limit  of  a  relatively  stiff  matrix  can  be  analysed  by  observing  that  aoz  while  /? 

^  that  Eq.  (3.2.9c)  in  the  limit  Exn  oo  therefore  yields  the  limiting  solution  tj  — >  (2cr)' 

(first  and  last  terms  of  Eq.  (3.2.9c)  dominant).  Thus  laprj  when  0. 

Figure  3b  shows  this  approximation  for  three  values  ofy0<  1.  While  the  limit  is  approached  only 
for  ap  so  small  that  lapr]  is  also  small,  it  can  nevertheless  be  used  as  the  basis  for  engineering 
estimates  of  the  condition  for  significant  inertial  effects,  i.e.,  the  first  significant  departure  of 
lapri  from  unity.  The  construction  of  Fig.  3b  suggests  that  inertial  effects  will  be  significant 

when  ap<  A{ap)u  where  {ap)x  is  the  value  of  the  product  a/?  at  which  the  equation 
=  1  is  satisfied;  i.e.. 


0^<1).  (3.2.17) 


This  leads  to  the  criterion  for  significant  inertial  effects  for  composites  with  relatively  stiff 
matrices  that 


<2| 


Iz/ 

/ 


n3 


ElE^  R 


T  C„ 


iP<\) 


(3.2.18) 


The  criteria  of  Eqs.  (3.2.16)  and  (3.2.18)  coincide  when  p=  1.  Since  laprf  <  1  always,  inertial 
effects  increase  the  stiffness,  dpidu,  of  the  traction  law  for  loads  that  increase  linearly  in  time. 

Conclusion: 


Some  analytical  results  have  been  presented  for  the  problem  of  bridging  by  the  mechanism  of 
fibre  pullout  when  the  inertia  of  the  fibre  and  the  matrix  are  taken  into  account.  With  this  model, 
we  also  calculated  simple  criteria  for  significant  inertial  effects  in  the  bridging  mechanism  in 
representative  mode  I  crack  propagation  problems.  We  find  that  the  inertial  effects  in  the 
bridging  mechanism  will  often  be  significant  for  a  matrix  crack  propagating  dynamically  in  the 
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Steady  state  ACK  limit  in  a  brittle  matrix  continuous  fibre  composite  or  delamination  in 
laminates  reinforced  by  through-thickness  stitching  or  rods. 

We  also  find  that  significant  inertial  effects  are  favored  by  low  fibre  volume  fi'action,  low 
fnction  stress,  low  matrix  bar  wave  speed,  and  low  fibre  modulus;  and  high  fibre  diameter  and 
high  matrix  modulus.  If  the  fibre  modulus  is  high  enough  (relative  to  the  matrix  modulus),  the 
criterion  for  inertial  effects  becomes  independent  of  fibre  modulus.  The  matrix  density  enters 
the  criterion  only  through  the  matrix  bar  wave  speed.  The  criterion  is  always  independent  of  the 
fibre  density. 

Finally,  we  also  find  that  for  pullout  or  bridging  stresses  that  rise  linearly  in  time,  the 
instantaneous  crack  displacement  is  less  in  the  presence  of  inertial  effects  than  it  would  be  under 
static  loading  to  the  same  bridging  stress.  However,  solutions  for  pullout  from  a  rigid  matrix 
suggest  that,  if  the  bridging  stress  rises  rapidly  and  is  then  held  at  a  constant  value,  the  crack 
displacement  when  all  particle  motion  finally  stops  will  be  greater  than  it  would  have  been  under 
static  loading  to  that  stress  level.  Thus  regimes  of  both  hardening  and  softening  of  the  bridging 
traction  law  due  to  inertial  effects  can  be  expected  in  bridged  crack  problems. 

We  also  formulated  fast  numerical  solutions  for  general  loading  cases  and  hence  have  a  scheme 
to  compute  the  bridging  law  for  fairly  general  conditions  [4].  This  will  be  required  to  solve 
large-scale  bridging,  dynamic  crack  problems  to  self-consistency  (bridging  tractions  unknown  a 
priori).  The  history  dependent  bridging  law  thus  computed  will  form  the  basic  ingredient  in  the 
computation  of  the  structural  response  of  the  composite  under  dynamic  loading  conditions. 
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3.4  Mechanics  of  Crack  Growth  in  Through-Thickness  Reinforced  DCB 
Specimens 

In  this  portion  of  the  research  [5,6],  the  role  of  bridging  by  stitches  or  rods  on  steady  state 
dynamic  crack  growth  is  analyzed.  In  particular,  the  dynamics  of  crack  growth  for  through 
thickness  reinforced  double  cantilever  beam  (DCB)  specimens  are  examined.  The  DCB 
specimen  is  widely  used  in  Mode  I  fracture  toughness  tests  of  polymers  and  composites. 
Certain  solution  characteristics  as  a  result  of  the  bridging  are  first  identified.  In  the  next  section, 
we  look  into  detail  into  the  energetics  of  crack  growth  for  a  through  thickness  reinforced  DCB 
specimen  loaded  by  a  flying  wedge.  The  double  cantilever  beam  (DCB)  specimen  loaded 
dynamically  by  a  flying  wedge  offers  a  relatively  simple  experimental  approach  to  analyzing  the 
mode  I  dynamic  delamination  problem.  The  test  is  especially  attractive  for  studying  the  bridging 
effects  supplied  by  through-thickness  reinforcement  (e.g.,  stitches  or  rods)  in  laminates.  The  role 
of  bridging  by  stitches  or  rods  on  steady  state  dynamic  crack  growth  was  computed  by  solving 
the  bridged  crack  problem  within  the  framework  of  beam  theory. 


Beam  Theory  Formulation  and  Solution  Characteristic.^' 

The  equations  of  motion  for  the  beam  element  are; 
aN_  d^u 

dK~  (3.4.1a) 


§-p(w,,)b=  pbh^ 

ax  Q^2 


ax 


(3.4.1b) 

(3.4.1c) 


where  u(x,t)  and  w(x,t)  are  the  in-plane  and  transverse  displacements  of  the  neutral  plane 
respectively,  ^(x,t)  is  the  clockwise  rotation  of  the  cross-section,  t  is  the  time  variable,  N  is  the  axial 
force,  Q  is  the  shear  force,  M  is  the  bending  moment,  2h  is  the  total  thickness  of  the  DCB  specimen, 
b  is  the  width  of  the  specimen,  p  is  the  density,  I  (=  bhVl2)is  the  moment  of  inertia  and  p(w,t)  is  t.he 
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bridging  traction  corresponding  to  the  opening  mode.  The  time  dependent  bridging  traction  p 
corresponding  to  the  opening  mode  is  assumed  to  depend  only  on  the  transverse  displacement  w. 


For  a  Timoshenko  beam,  the  equations  of  motion  for  steady  state  cracking  reduce  to; 


gV  _ cf _ 12  R 

ax'*  (R-cf)(l-cf)  h^  ^ 

1  1  a^p(w,X)  1  12R  , 

- 5 - ; - 5 — ~  + - 5 - 0 - rP(w,X)  =  0 

(R-cf)Eh  aX^  (R-cf)(l-c,^)  Eh^ 

_  p{w,X)  (R-cj)  d^w 

^  ~  REh  'r 


(3.4.2a) 


(3.4.2b) 


(3.4.2c) 


where  X  =  x-  vt,  c]  ■=  pv^ !  E  and  R=kGI  E  and  where  v  is  the  (constant)  velocity  of  the 
crack  tip.  In  the  limit  of  steady  state  crack  velocity  v— >0,  we  obtain  the  equations 
corresponding  to  the  static  case. 


Similarly,  it  is  easy  to  show  that  for  and  Euler-Bernoulli  beam  (where  both  the  shear  deformation 
and  rotational  inertia  is  ignored),  equations  of  motion  for  steady  state  cracking  reduce  to; 


(3.4.3a) 


12  cf  aV 

h^  dX} 


+ 


12 

Eh^ 


p(w,X) 


=  0 


(3.4.3b) 


where  X  -  x  -  v  t,  v  is  the  (constant)  velocity  of  the  crack  tip,  c\  =  pv^ !  E  and  R=  kGIE  . 
We  shall  now  consider  a  linear  bridging  law  of  the  following  type; 

p(w,X)  =  p„  + J3^w  (3.4.4). 

A  linear  bridging  law  of  the  type  above  is  a  simple  but  realistic  constitutive  law.  In  the  results 
that  follow,  we  non-dimensionalize  the  length  variables  by  the  laminate  thickness  h.  Therefore, 
w  shW ^  u  =  hU  and  X  =  h  ^ .  We  now  analyze  the  solution  characteristics. 
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For  a  linear  bridging  law  in  Equation  3.4.4,  t/(^)=0.  The  transverse  displacement  obeys: 


'  ^  ^  +  b^W  +  d'^  =0 


Also,  for  the  Timoshenko  beam  the  clockwise  rotation  (|)  is: 

^  _  {R-cl)  d^W 

a^  RE  R 

The  coefficients  3,  b  and  d  in  Equation  3.4.5a  are: 


(3.4.5a) 


(3.4.5b) 


For  Timoshenko  beam. 

P-1 


12  cf  R 


(l-c;)(i?-c/)  {R-c^)E 


b  = 


d  = 


R 


np^h 


\\-c^){R-c^)  E 


R 


12/^0 


'(l-c/)(i?-cf)  E 


(3.4.6a) 


For  Evler-Bernoiilli  (E-B)  beam: 


(3.4.6b) 


and  where  c]  -  pv^ ! E  ^  c]  =  pv^  I{k G)  and  R  =  kG  /E.  In  the  limit  that  the  crack  velocity 
approaches  zero  (  v  =  c,  =  — >  0),  we  retrieve  the  static  case  of  Roberta  and  Cox  [17,18]. 


The  general  solution  to  equation  3.4.5a  is: 


d^  - 

m)  =  ~  +  K,e  V  ^ 

b 

- 

+  K,e  2 


"  +K,e 


V  2  2^ 


Ab^  4 


(3.4.7) 
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There  are  three  regimes  to  the  solution  behavior  and  these  are  identified  below.  In  the  equations 
below,  we  have  introduced  the  dimensionless  quantities; 


•  Case  1:  3^  <  0  and  3“*  >  4b^  =>  Exponential  behavior 

-  For  Timoshenko  beam,  this  is  true  provided: 

and  3(S(l-c,=  )-c’)=>S(l-c;)(/!-c?)  (3A8a) 

E  1  +  0 

-  For  the  E-B  beam,  the  above  condition  is  never  satisified. 

•  Case  2;  3“  ^  0  and  3"*  >  4b^  =>  Oscillatory  and  non-decaying  behavior 

-  For  Timoshenko  beam,  this  is  true  provided: 

and  2{S{\-c^)-c^y>S{\-c^){R-c^)  (3.4.8b.l) 

E  1+0 

-  For  the  E-B  beam,  this  condition  is  satisified  when: 

2 

>  2^^  (3.4.8b.2) 

E 

•  Case  3 :  3"*  ^  4b^  =>  Oscillatory  with  exponential  decay  behavior 

-  For  the  Timoshenko  beam,  this  is  true  when; 

3{S(l-cf)-cfy  <S(l-cfXR-cf)  (3.4.8C.1) 

-  For  the  E-B  beam,  this  is  true  when; 

2yf^  (3.4.8C.2) 

E 

The  arbitrary  constants  Ki,  K2,  K3,  and  K4  in  Equation  15  are  defined  through  the  associated 
boundary  and  continuity  conditions.  The  conditions  determined  above  give  us  insight  into  when 
dynamic  effects  can  significantly  alter  the  mechanisms  of  deformation  and  the  resultant  bridging 
phenomena.  For  instance,  if  the  crack  tip  velocity  exceeds  the  conditions  prescribed  in  either 
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(3.4.8a),  oscillatory  displacement  fields  will  be  introduced  in  the  wake  of  the  crack.  These 

multiple  oscillations  are  of  a  very  different  form  to  the  crack  face  interpenetration  caused 

sometimes  by  bridging  effects  in  the  static  case.  When  such  oscillations  result  from  dynamic 

effects,  the  mechanics  of  bridging  and  the  efficacy  of  through  thickness  bridging  ligaments  on 

the  energetics  of  crack  growth  will  be  considerably  altered.  For  example,  stick-slip  propagation 

modes  become  likely. 


Wedge-Loaded  Double  Cantilever  Beam 


Figure  4:  Schematic  of  through-thickness  reinforced  DCB  specimen  loaded  with  a  flying  wedge 

We  now  concentrate  on  the  mechanics  of  delamination  for  the  Euler  Bernoulli  beam, 
where  both  shear  deformation  and  rotational  inertia  can  be  ignored.  Therefore  the  constants  p,  b 
and  d  referred  to  hereafter  are  the  ones  presented  in  Eqn.  3.4.6b.  The  double  cantilever  beam 
(DCB)  specimen  loaded  dynamically  by  a  flying  wedge,  of  constant  velocity,  offers  a  relatively 
simple  experimental  approach  to  studying  the  mode  I  dynamic  delamination  problem  (Fig.  4). 
The  test  is  especially  attractive  for  studying  the  bridging  effects  supplied  by  through-thickness 
reinforcement  (e.g.,  stitches  or  rods)  in  laminates.  In  figure  4,  2a  is  the  wedge  angle,  I  is  the 
distance  between  the  wedge  and  the  crack  tip  and  ao  is  the  length  of  the  bridging  zone.  In  non- 

dimensional  form,  I  =  h  L,  and  =h  A^.  Due  to  symmetry,  we  focus  our  attention  only  on  the 
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top  half  of  the  DCB  specimen.  A  detailed  analysis  of  the  role  of  the  bridging  on  dynamic 
delamination  behavior  is  presented  in  this  section.  The  analysis  is  valid  when  the  bridging 
constituent  has  a  linear  bridging  relationship  (as  presented  in  Eqn.  3.4.4). 


For  the  unbridged  portion,  the  deflection  profile  (w„  =hW^)  is  obtained  by  setting  b  =  d  =  0. 
Therefore: 


dX 

d^* 

o 

11 

for  {-L<^<-A,) 

(3.4.9a) 

d^W 

— ; — V  b'W  +  d"  =0 

for  {-A^  <^<0) 

(3.4.9b) 

The  relevant  boundary  conditions  are: 

^(#  =  0)  =  0,  W'(4  =  0)  =  0, 

W:X^  =  -L)  =  -a,  w:(^  =  -L)  =  0  .  (3.4.10) 

The  deflection  profiles  should  satisfy  the  continuity  conditions  at  the  end  of  the  bridging  zone 
(^  =  -Aq).  The  continuity  conditions  are: 


Wi^  =  -A,)  =  W„(4  =  -A,),  =  -A,)  =  =  -A,), 

.  (3.4.11) 

The  governing  equation  (3.4.9)  together  with  the  boundary  conditions  (3.4.10)  and  the  continuity 
conditions  (3.4.11)  will  completely  determine  the  deflection  profile  of  the  beam.  The  bridging 
zone  length  (Ao)  will  be  dictated  by  the  critical  crack  opening  displacement  (vi'^  =hW^)  required 
for  failure  of  the  bridging  ligament.  Therefore,  Ao  is  obtained  by  solving: 

^(4  =  -Ao)  =  fV,  .  (3.4.12) 


The  energy  release  rate  (  ),  determined  through  the  total  energy  balance  is: 


'Total 


_  ifaj„  au.  auj 


B 


da 


da  da 


(3.4.13) 


where  Uext  is  the  work  done  by  the  applied  load,  Us  is  the  strain  energy,  Ur  is  the  kinetic  energy, 
B  is  the  uniform  width  of  the  DCB  specimen,  and  a  is  the  crack  length.  For  steady  state  crack 
extension  a  =  v  t,  where  v  is  the  crack  velocity  and  t  is  time.  Further,  under  steady  stair 
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delamination,  there  is  no  net  change  in  strain  energy  (Us  =  0)  upon  crack  extension  since  both 

the  material  ahead  of  the  crack  tip  and  the  material  in  the  wake  of  the  wedge  have  zero  strain 
energy.  However,  for  small  a: 

o'Ui.  .5^1  2'^  IT  T, 

"aT  J  ^  Sin^(a)«  Bphv^a^  (3.4.14a) 

au  , 

■  =  2  .Applied  Force,  v 

at 


=  -2Sin(a)v 


Bv^^ 

6  dC 


f,  =  -L 


(3.4.14b) 


where  m  is  the  mass  in  the  wake  behind  the  wedge,  v„  is  the  normal  velocity  and  Q  is  the  shear 
force  in  the  beam.  Therefore, 


^Total 


aEh  a^W, 
6 


phv^a^ 


(3.4.15) 


In  addition,  by  application  of  the  dynamic  J-integral,  the  energy  released  at  the  crack  tip  is 
related  to  the  bending  moment  M  by  : 


12 

Eh 


[m 

L  _  Eh 

■N 

L  4=0  . 

^  "  12 

(3.4.16) 


In  general,  due  to  crack  bridging,  AG  -  -  G^-jp  ^  0.  The  energy  difference,  AG  ,  is  the 


extra  work  done  in  fracturing  the  bridging  ligaments.  In  the  absence  of  crack  bridging,  AG=  0 
^Totai  =  =  Go  is  given  by: 


G  =  P  ^ 

°  Tan^(X/2) 


(3.4.17) 


With  this  formalism,  detailed  calculations  of  the  deflection  profile  and  crack  energy  release  rate 
was  computed  for  the  linear  bridging  law  [5].  Analytical  results  were  obtained  for  steady  state 
crack  propagation  conditions.  The  key  results  are  summarized  below: 

•  The  magnitude  of  the  bridging  traction,  which  is  controlled  by  the  properties  of  the  through 
thickness  reinforcement,  plays  a  significant  role  in  the  crack  growth  behavior.  For  steady 
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State  crack  growth  conditions,  different  regimes  of  the  solution  behavior  have  been  identified 
which  would  correspond  to  different  crack  propagation  characteristics.  Critical  velocities 
that  demarcate  the  various  solution  regimes  have  been  identified. 

•  Steady  state  crack  growth  is  attainable  for  this  loading  configuration  provided  certain 
conditions  are  satisfied.  Regions  of  stable  crack  growth  as  a  function  of  the  material 
properties  of  the  through  thickness  reinforcement,  the  size  of  the  DCB  specimen  and  the 
velocity  of  the  wedge  have  been  identified.  This  provides  guidelines  for  design  of 
experiments  to  probe  the  efficacy  of  bridging  on  improving  the  dynamic  fracture  toughness 
of  through  thickness  reinforced  structures. 


3.5  Dynamic  Crack  Energy  Release  Rate  for  Orthotropic  Laminates 

In  this  portion  of  the  work,  we  compute  explicit  relations  for  the  dynamic  crack  energy  release 
rate  and  crack  tip  stress  intensity  factors  in  a  form  most  amenable  to  computational  work  for  a 
crack  propagating  along  any  principal  axis  in  an  orthotropic  material  [6],  The  results  presented 
in  this  form  may  be  particularly  useful  for  work  on  delamination  cracks  in  laminated  composites, 
which  are  orthotropic  in  many  important  applications.  Although  Yang  and  co-workers  [13] 
presented  results  for  a  mixed  mode  crack  in  an  orthotropic  material,  it  is  in  a  form  that  is  not  best 
suited  for  immediate  computation.  By  an  adaptation  of  the  Stroh  method  of  analysis  used  by  Wu 
[14],  the  results  of  Yang  et  al.  [13]  are  rendered  here  in  a  different  form  better  suited  to 
numerical  evaluation.  The  role  of  orthotropic  anisotropy  in  determining  the  velocity  dependence 
of  Gd  is  then  illustrated  for  delamination  cracks  in  typical  ceramic,  metal  and  polymer  matrix 
laminated  composites. 

The  relation  between  the  energy  release  rate  and  the  crack  tip  stress  intensity  factor  is  central  to 
fracture  analysis.  It  appears,  for  example,  in  the  derivation  of  stress  intensity  factors  from 
applications  of  the  J-integral  and  in  the  derivation  via  energy  arguments  of  integral  equations  for 
bridged  cracks  from  weight  functions.  Freund  [15]  showed  that  the  dynamic  crack  energy 
release  rate  for  an  extending  crack  in  an  elastic  body  can  be  written  as  the  modified  J-integral, 
which  leads  to  the  G-k  relationship  in  the  following,  convenient  form; 
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Gd  =  L-‘k(t) 

where  Gd  is  the  dynamic  crack  energy  release  rate  and  k(t)  is  the  instantaneous  crack  tip  stress 
intensity  factors  (see  Eq.  4.19  in  Wu  [14]  and  references  therein).  The  elements  of  the  U*  matrix 
are  universal  functions,  in  the  sense  that  they  are  independent  of  the  details  of  the  applied 
loading  or  on  the  configuration  of  the  body  being  analyzed.  We  calculated  the  L'^  matrix  when 
the  material  possesses  orthotropic  symmetry  and  when  the  crack  is  extending  along  one  of  the 
principal  axis  [6].  In  this  case,  the  off-diagonal  elements  of  the  matrix  are  zero  and  the  diagonal 
elements  of  the  L'*  matrix  is  just  dependent  on  the  elastic  properties  of  the  orthotropic  medium 
and  the  instantaneous  value  of  the  crack  tip  velocity.  For  an  orthotropic  material  with  the  crack 
propagating  along  the  1  direction  (and  direction  3  being  the  plane  strain  direction),  elements  of 
the  L'*  matrix  are: 


^-(r5+ri/4) 

(rA-rsrs) 


0 

-(rs-ri/A) 

(rl-rsrs) 


0  0 


0 

0 

>'6; 


where  the  y’sare  presented  in  detail  in  the  appendix  paper  [6],  We  can  also  define  the 
generalized  Rayleigh  wave  function  R(v)  as 


^(v)  =  rhrs/s 


and  the  Rayleigh  wave  speed  (Vr)  is  obtained  by  setting  R(vr)  =  0. 

The  diagonal  elements  L  and  L  '2,2  were  computed  numerically  for  cracks  propagating  along 
different  principal  axes  in  various  representative  composites.  The  L**  matrix  elements  are 
functions  of  the  elastic  properties  and  the  instantaneous  crack  velocity.  Each  diagonal  element 
approaches  the  corresponding  static  value  as  v  — >  0  and  has  the  property  of  0\j(\\  —  v)~*  ]  as 
’  where  Vr  is  the  Rayleigh  wave  speed.  In  addition,  we  find  that  the  that  the  variation  in 

the  normalized  crack  energy  release  rate  (L‘*  /  [L**]static)  as  a  function  of  the  normalized  crack 
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velocity  v/vr  can  be  rationalized  in  terms  of  the  orthotropy  parameters  A,  and  p  (Suo  et 
al.,[16]) .  For  the  same  normalized  crack  speed,  the  normalized  mode  I  contribution  /  [L* 
\,i]static)  to  the  dynamic  crack  energy  release  rate  increases  as  X  monotonically  increases  and  as 
p  monotonically  decreases.  For  the  same  normalized  crack  speed,  the  normalized  mode  II 
contribution  (L‘S,2/  [L**2,2]static)  to  the  dynamic  crack  energy  release  rate  increases  as  X  increases 
and  as  p  decreases.  However,  as  opposed  to  the  mode  I  case,  the  variation  on  X  and  p  is  much 
smaller.  In  addition,  we  observe  that  although  the  energy  release  rate  shows  a  monotonic 
dependence  on  X,  the  dependence  on  p,  for  small  p,  is  not  monotonic. 

To  conclude,  we  have  provided  an  analytical  form  for  easy  estimation  of  the  dynamic  crack  energy 
release  rate,  in  terms  of  the  crack  tip  stress  intensity  factors,  has  been  presented  for  cracks 
propagating  along  any  of  the  principal  axis  in  orthotropic  material  systems.  The  dynamic  crack 
energy  release  rate  depends  on  the  magnitude  of  the  two  orthotropic  parameters  and  the 
instantaneous  crack  tip  velocity.  The  dynamic  crack  energy  release  rate  is  0[(v^  -  v)”’  ]  as  v  — >  v,. , 

and  where  Vr  is  the  Rayleigh  wave  speed.  The  variation  in  the  dynamic  crack  energy  relase  rate  for 
orthotropic  materials  can  be  rationalized  in  terms  of  the  two  orthotropy  parameters  X  and  p. 
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5.  COLLABORATIONS 

We  continue  to  collaborate  with  Dr.  Irene  Beyerlein  of  Los  Alamos  National  Labs  and  Dr. 
Roberta  Massabo  of  University  of  Genova,  Italy  on  modeling  the  mechanisms  of  deformation 
under  dynamic  loads. 

Prof  Ares  Rosakis  of  Caltech,  with  a  sub-contract  under  this  program  (~$20K),  is  conducting 
experiments  for  detailed  study  of  the  dynamic  fiber  pull-out  /  push-in  process  to  provide  data 
with  which  the  validity  of  the  pull-out  model  will  be  assessed.  Through  the  use  of  a  carefully 
designed  2D  model  material  system  and  high-speed  diagnostics,  a  complete  quantitative 

description  of  process  can  be  obtained  which  makes  direct  comparison  with  the  theoretical  model 
possible. 
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ABSTRACT 

Incrti&l  effects  in  the  mechflnism  of  fibre  pullout  during  the  dynsniic  propagation  of  a 
bridged  crack  are  examined  by  reposing  simple  shear  lag  models  of  pullout  as  problems 
of  dynamic  wave  propagation.  The  only  coupling  considered  between  the  fibres  and  the 
matrix  is  uniform,  rate  independent  friction  —  no  debond  energy  is  included.  Analytical 
solutions  are  found  to  the  problem  of  the  coupled  waves  in  the  fibres  and  the  matrix  that 
propagate  away  from  the  fracture  plane  of  the  bridged  crack  as  the  bridging  tractions 
increase  with  time.  These  solutions  yield  the  time-dependent  relationship  between  the 
crack  opening  displacement  and  the  bridging  traction.  Engineering  criteria  for  inertial 
effects  being  significant  are  deduced  by  comparing  the  dynamic  bridging  traction  law 
with  its  counterpart  for  static  loading,  which  is  recovered  as  a  limit  of  the  dynamic  case. 
The  criteria  are  evaluated  for  two  crack  cases:  the  asymptotic  limit  of  a  long,  fully 
bridged  matrix  crack  propagating  unstably  through  a  fibrous  composite  under  remote 
tension;  and  a  delamination  crack  bridged  by  stitches  or  rods  that  is  loaded  by  a  flying 
wedge  splitting  a  double  cantilever  beam.  In  both  cases,  the  rate  of  increase  of  the  crack 
opening  displacement  appears  to  be  sufficient  for  inertial  effects  to  be  pronounced  in  the 
bridging  (pullout)  mechanism.  Expected  trends  of  the  significance  of  inertial  effects  with 
material  and  geometrical  parameters  are  identified. 
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1.  Introduction 

The  pullout  mechanism  is  the  fundamental  source  of  toughening  and  fracture  resistance 
in  many  composites.  In  brittle  matrix  composites,  reinforcing  fibres  that  are  weakly 
bonded  to  the  matrix  can  survive  the  passage  of  matrix  cracks,  across  which  they  then 
provide  bridging  tractions  that  shield  the  crack  tip  from  the  applied  load.  The  main 
mechanism  of  load  transfer  from  the  fibres  to  the  matrix  is  interfacial  friction.  In 
polymeric  laminates  reinforced  through  the  thickness  by  stitches  or  rods,  analogous 
pullout  phenomena  are  observed,  but  on  the  scale  of  fibre  tows,  which  may  be  1  mm  in 
diameter,  rather  than  on  the  scale  of  individual  fibres  (10  -  100  pm).  Scale 
considerations  aside,  the  mechanics  of  pullout  are  very  similar  in  the  two  cases.  Other 
systems  in  which  bridging  entities  are  coupled  to  a  matrix  by  friction  include  so-called 
self-reinforced  polycrystalline  ceramics,  in  which  elongated  grains  bridge  cracks;  and 
ceramic  layered  systems,  in  which  fractured  layers  slide  past  one  another  during  failure. 

For  mode  I  cracks,  the  shielding  effect  created  by  the  pullout  phenomenon  can  be 
summarised  by  a  bridging  traction  law  that  relates  the  stress  in  the  bridging  entities  at  the 
fracture  plane,  T,  to  the  crack  opening  displacement,  2ti.  To  a  very  good  approximation 
in  many  cases,  the  tractions  in  the  bridging  entities  can  be  replaced  by  an  equivalent 
continuous  traction,  p,  that  is  to  be  applied  to  the  entire  bridged  interval  of  the  fracture 
surfaces  (e.g.,  Cox  and  Marshall,  1994).  Predicting  the  traction  law,  p{ti),  becomes  one 
of  the  central  problems  of  crack  bridging  theory. 

The  mechanics  of  pullout  and  the  resulting  traction  law  have  been  much  studied  and  are 
well  understood  for  static  loading.  Simple  anal34;ical  forms  are  available  for  p(ti)  when 
the  frictional  coupling  of  the  reinforcement  to  the  matrix  is  uniform  and  slip  extends  over 
distances  that  are  large  compared  to  the  reinforcement  diameter  (Marshall,  Cox,  and 
Evans,  1985;  McCartney,  1987).  In  this  limit,  which  is  a  common  case  in  ceramic 
composites  and  textile  polymeric  composites,  the  shear  lag  model  of  load  transfer 
between  the  reinforcement  and  the  matrix  is  accurate.  Simple  extensions  of  models  of 
this  class  are  also  available  to  deal  with  small  but  nonzero  levels  of  the  work  required  to 
debond  the  reinforcement  from  the  matrix  prior  to  slip  (Hutchinson  and  Jensen,  1990). 

Given  the  relationship,  p(iu),  the  characteristics  of  crack  propagation  can  be  computed  by 
solving  a  bridged  crack  problem.  In  composites  in  which  the  bridging  mechanism  is 
most  effective,  the  zone  of  bridging  can  be  comparable  to  the  crack  length  and  much 
larger  than  features  such  as  notches  (Bao  and  Suo,  1992;  Cox  and  Marshall,  1994)  or,  in 
the  case  of  delamination  cracks,  the  laminate  thickness  (Jain  and  Mai,  1995;  Massabo  and 
Cox,  1999).  Crack  propagation  then  does  not  follow  Linear  Elastic  Fracture  Mechanics, 
in  the  sense  that  there  is  no  single  material  parameter  such  as  toughness  that  correlates 
with  the  crack  growth.  Instead,  the  bridged  crack  problem  is  one  of  large  scale  bridging 
and  the  characteristics  of  propagation  show  features  peculiar  to  the  form  of  the  traction 
law, /?(//)  [Cox,  1991]. 

This  paper  extends  existing  models  of  the  mechanics  of  pullout  to  high  loading  rates.  An 
approach  to  evaluating  a  traction  law  that  takes  account  of  the  inertia  of  the  reinforcement 
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and  the  matrix  is  formulated  as  a  direct  extension  of  the  elementary  static  loading  model 
of  McCartney  (1987)  and  Marshall,  Cox,  and  Evans  (1985).  The  chosen  base  model  for 
static  loading  has  proven  consistent  with  experiments  in  many  material  systems  and  has 
been  the  foundation  of  major  advances  in  understanding  damage  in  composites.  The 
spirit  of  the  present  work  is  to  seek  equivalent  insight  into  dynamic  damage  by 
incorporating  the  influence  of  inertia  into  the  simplest  credible  model.  Thus  not  all 
aspects  of  the  micromechanics  of  pullout  that  could  be  important  in  some  cases  will  be 
addressed.  Instead,  attention  will  be  focused  on  identifying  a  characteristic  time  for  the 
frictional  pullout  problem  that  will  allow  rapid  assessment  of  when  inertial  effects  will  be 
important. 

2.  Idealisation  of  the  Bridged  Crack  problem 

The  crack  propagation  and  pullout  problems  are  depicted  schematically  in  Fig.  la.  A 
matrix  crack  propagates  on  the  plane  z  =  0  and  is  bridged  by  intact  fibres  in  its  wake 
(For  simplicity  of  expression,  the  term  fibre  from  here  on  will  be  used  to  refer  to  bridging 
entities  of  any  kind,  including  stitches,  rods,  and  bridging  grains.)  Upon  the  passing  of 
the  crack  tip,  a  debond  crack  propagates  along  the  length  of  each  fibre  away  from  the 
fracture  plane.  Propagation  of  the  debond  crack  is  governed  by  the  fracture  energy 
associated  with  the  separation  of  the  matrix  and  the  fibre  at  the  debond  crack  tip  and  the 
work  done  against  friction  in  displacing  the  debonded  fibre  along  its  axis.  In  many 
composites,  the  debond  energy  is  small  and  pullout  is  dominated  by  friction  over  much  of 
the  range  of  pullout  displacements.  Therefore,  in  this  first  model  the  debond  energy  will 
be  ignored. 


For  static  loading,  the  crack  propagation  problem  can  be  idealised  by  replacing  the 
process  zone  by  elastic  composite  material  down  to  the  fracture  plane  and  representing 
the  phenomena  within  the  process  zone  by  bridging  tractions  applied  continuously  on  the 
fracture  surfaces  (Fig.  lb).  The  bridging  tractions,  p,  are  related  to  the  axial  stress,  T,  in 
the  fibre^  at  the  fracture  plane  by/;  =/r,  where /is  the  area  fraction  of  the  fibres  on  the 
plane  z  -  0.  For  aligned  continuous  reinforcement,  /  is  also  the  volume  fraction  of  the 
i.  crack  opening  displacement,  2//,  in  the  idealisation  should  be  defined  as 

the  difference  in  the  actual  displacement  evaluated  across  the  process  zone  and  the 
isp  acement  that  would  be  expected  if  the  material  in  the  process  zone  was  elastic. 

The  traction  law, /?(m),  is  derived  by  considering  the  micromechanics  of  the  phenomena 
occumng  within  the  process  zone,  which  is  to  say  the  micromechanics  of  frictional 
s  1  in^  he  micromechanical  problem  can  be  represented  by  a  small  volume  of  material, 
e.g.,  the  material  bounded  by  one  of  the  dotted  rectangles  in  Fig.  la.  The  traction 
boundary  conditions  for  the  representative  volume  are  as  follows.  At  z  =  0,  the  matrix  is 
traction-free,  while  the  fibres  sustain  the  axial  traction  T.  At  z  =  4,  the  strain  in  the  fibres 
an  t  e  matrix  must  equal  the  average  strain  in  the  composite  adjacent  to  the  process 
zone  (z  >  4).  Shear  tractions  may  arise  along  the  vertical  boundaries  of  the  representative 
volume  (parallel  to  z),  but  these  are  neglected. 
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In  the  problem  of  a  crack  propagating  under  static  loading,  the  traction  at  z  =  4  for  the 
representative  volume  associated  with  a  fixed  volume  of  material  is  assumed  to  rise  from 
zero  when  the  material  is  immediately  behind  the  crack  tip  to  increasing  values  as  the 
material  passes  further  back  into  the  wake  (Marshall  et  al,  1985).  The  tractions  are  not 
truly  zero  right  at  the  crack  tip,  since  the  fibres  are  not  stress-free,  but  assuming  they  are 
leads  to  reasonable  results  for  the  bridged  crack  problem  as  long  as  the  fibres  remain 
intact  over  sufficiently  long  distances  into  the  crack  wake.  Then  shielding  of  the  crack 
tip  is  dominated  by  the  bridging  tractions  acting  in  the  further  crack  wake,  where  the 
boundary  conditions  are  correct. 

Further  details  of  the  analysis  of  the  micromechanical  problem  may  be  found  in 
McCartney  (1987)  and  the  appendix  to  Marshall  et  al.  (1985). 

In  dynamic  loading,  the  boundary  conditions  at  the  boundary  of  the  process  zone  involve 
displacement  and  displacement  rates  as  well  as  stress  or  strain  conditions. 

3.  The  Micromechanics  of  Dynamic  Pullout 

When  friction  is  the  only  active  force  of  resistance,  the  static  pullout  problem  reduces  to 
an  idealization  in  which  fibers  are  pulled  out  of  a  half-space  to  which  they  are  not 
initially  bonded.  The  dynamic  problem  is  one  of  wave  propagation  along  a  fibre  in  the 
presence  of  frictional  retardation  and  with  no  debond  crack  tip  or  crack  tip  field  to  be 
considered. 

The  archetypal  problem  is  illustrated  in  Fig.  2.  A  representative  volume  consists  of 
cylindrical  fibres  of  radius  R  and  volume  fraction  /  embedded  in  a  matrix  (z  >  0).  The 
fibre  and  the  matrix  have  axial  Young’s  moduli  E{  and  Em  and  densities  p[  and 
respectively.  The  axial  displacement,  strain,  and  stress  of  the  fibre  and  the  matrix  are 
denoted  ttt  and  Um,  fr  and  em,  and  Of  and  am,  respectively.  The  axial  displacements  will  be 
assumed  to  be  the  only  nonzero  displacement  components  induced  by  loading  and  to  be 
uniform  across  any  section  of  the  fibre  or  the  matrix.  Thus  the  displacement,  strain,  and 
stress  in  the  fibre  and  the  matrix  are  functions  of  z  and  t  only.  These  are  the  usual 
assumptions  of  shear  lag  theory  with  the  simplest  conditions  of  elasticity  (Poisson’s 
effect  omitted).  They  are  consistent  with  assuming  that  the  friction  forces  are  constant 
(unaffected  by  fibre  contraction  or  dilation  due  to  axial  stresses). 

There  is  no  initial  bond  between  the  fibres  and  the  matrix.  The  fibres  are  coupled  to  the 
matrix  by  friction  tractions,  Tf,  which  are  assigned  the  following  properties.  If  relative 
motion  exists  between  the  fibres  and  the  matrix,  then 


rf=  r 

("f  <"m) 

(la) 

Tf  =  -r 

(»f  >"m) 

(lb) 
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where  a  dot  indicates  time  differentiation,  r  is  a  positive  constant  and  rr  >  0  indicates 
friction  tractions  acting  on  the  fibres  in  the  positive  z  direction.  When  the  fibres  and  the 
matrix  are  not  in  relative  motion,  the  friction  tractions  may  support  stress  gradients  in  the 
fibres  and  the  matrix,  provided  that  the  required  magnitude  of  Tf  does  not  exceed  r.  Thus, 
by  equilibrium  considerations. 


dcTf  diir 


& 

R 


d<j„ 


dz 


-Pn 


diL 


dt 


2/ 


\-fR 


("f=Wm)-  (IC) 


Here  the  possibility  is  included  that  the  matrix  and  the  fibres  have  the  same  non-zero 
velocity  and  may  also  be  accelerating  together,  although  such  general  solutions  will  not 
be  exhibited  in  this  paper.  With  such  a  friction  law,  the  dynamic  wave  equations 
describing  those  parts  of  the  fibres  and  the  matrix  that  are  in  relative  motion  may  be 
written  approximately  as 


_  _  29x  1 


(2a) 


and 


A,  ,2/  $t  \  dhi„ 

1-/RE^  d  dt^ 

where  1  if  Wj.  and  0=  -\  otherwise;  and  Cf  and  Cm  are  the  bar  wave  velocities  in 
the  fibres  and  the  matrix,  given  by 


The  approximation  of  using  the  bar  wave  velocities  rather  than  the  longitudinal  wave 
velocities  in  the  z  direction  in  the  wave  equations  is  consistent  with  the  simplified 
treatment  of  stresses  and  strains  and  the  assumed  uniformity  of  the  friction  stress. 

Boundary  conditions  in  the  dynamic  case  are  as  follows.  At  the  fracture  plane, 

«f(0,  r)  =  0  &  CTmCO, /)  =  0  (z  =  0).  (3) 

At  the  boundary  of  the  process  zone  (limit  of  relative  fibre/matrix  motion), 

Wf  Wm  =  Wc  (z  =  4)  (4) 
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where  Uc  is  the  displacement  of  the  adjacent  intact  composite;  and  conditions  also  exist 
on  stress  or  strain  and  particle  velocities.  These  further  conditions  depend  on  the  nature 
of  the  loading  history,  which  can  be  expressed  as  the  function  £{t),  where  s  is  the  strain  in 
the  z  direction  in  the  intact  composite  adjacent  to  the  process  zone  boundary.  The 


bridging  traction,  p,  is  related  to  ^by 

p-  eE 

(5a) 

where  the  composite  modulus,  E,  is  given  by 

E  =fE{  +  (1  -f)Em 

(5b) 

In  the  depiction  of  Fig.  la,  the  process  zone  boundary  will  propagate  away  from  the 
fracture  plane  as  £{t)  rises.  In  the  case  to  be  considered  in  this  paper,  £(/)  will  be  assumed 
to  rise  continuously  and  monotonically  from  zero.  Placing  the  origin  of  time,  /  =  0,  at  the 
onset  of  nonzero  £(0.  the  location  of  the  zone  boundary  at  time  t  may  then  be  written^ 

4  =  7](t)Cmt  (6a) 

where  the  function  ri{t)  depends  on  material  and  geometrical  parameters  and  the  form  of 
s{t)  and  will  be  shown  to  be  bounded  by 

0  <  7(0  <  1  .  (6b) 

For  the  loading  conditions  considered,  the  additional  boundary  conditions  at  z  =  4  are; 

"f=«m=«c  (-  =  4)  (7a) 

and 

ff=£;n  =  £(0  iz  =  h).  (7b) 

If  the  load  history,  £(0,  possesses  discontinuities,  e.g.,  a  step  load,  then  discontinuities  in 
stress  and  particle  velocity  will  also  propagate  at  the  boundary  of  the  process  zone  (e.g., 
Achenbach,  1973;  see  also  Appendix  A). 

4.  Composite  Stress  Rising  Linearly  in  Time 

A  case  of  representative  interest  for  dynamic  bridged  crack  problems  and  for  which 
analytical  results  can  be  found  is  that  of  a  load  or  bridging  traction  that  increases  linearly 
in  time.  A  linearly  increasing  load  might  give  insight,  for  example,  into  bridging  effects 
in  a  specimen  in  which  substantial  bending  arises,  such  as  a  standard  double  cantilever 


'  It  is  easy  to  show  that  all  the  ensuing  results  for  linearly  rising  loading,  e  =  kt,  are  unchanged  if  7  is 
defined  instead  by  writing  4  =  and  thus  do  not  depend  on  the  relative  magnitudes  of  Cf  and  c^. 
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beam  delamination  specimen.  In  such  specimens,  the  crack  profile  is  often 
approximately  linear  and  the  rate  of  increase  of  the  bridging  tractions  at  any  point  might 
also  be  approximately  linear  if  the  crack  propagates  at  approximately  constant  speed. 
The  bridging  traction  at  a  particular  material  point  might  be  expected  to  rise  from  zero  as 
the  delamination  crack  first  passes  until  a  peak  value  is  reached,  perhaps  corresponding  to 
bridging  fibre  rupture.  Analytical  results  can  be  found  for  linearly  increasing  loads. 

Let 


sit)  =  kt 


(8) 


where  k  is  constant  and  all  displacements  and  boundary  tractions  are  zero  for  /  <  0.  In 
this  case,  0=  1  and  the  wave  equations  have  the  solutions 


k  2 

Wf  = - ca^  + 


■  [l  +  2r]a\ 


where  t]  is  independent  of  time  and  satisfies 
arf  +  Virf  +  Parj  - 14  =  0 

with  the  dimensionless  parameters  a  and  given  by 


(9a) 

(9b) 


(9c) 


TC 

a  =  ■ — 


& 


EfRk  '  i-fE^ 

Analysis  shows  that  Eq.  (9c)  has  only  one  real  root,  which  always  lies  in  (0,1). 
The  particle  velocities  and  the  strains  for  0  <  z  <  4  are  given  by 
Uf  =  [l  +  2Tja]/cz 
u^=[l  +  2rja]kl. 


(9d) 


Sf  = 


^  M 

l  +  2T]a 

<  4 

_ 

(10a) 

(10b) 

(10c) 


4 


(lOd) 
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The  matrix  particle  velocity  is  uniform  and  increases  in  proportion  to  4.  The  fibre 
velocity  rises  linearly  through  the  process  zone.  A  strain  concentration,  sds,  propagates 
in  the  fibre  behind  the  process  zone  front,  while  the  matrix  strain  in  the  process  zone  is 
always  less  than  s.  The  strain  distributions  and  velocities  beyond  the  process  zone,  z  >  4, 
need  not  be  specified,  provided  the  composite  strain  s  obeys  the  condition  Eq.  (7)  on  the 
process  zone  boundary. 

The  displacement,  ii,  to  be  used  in  defining  the  traction  law, /?(?/),  is  given  by  the  common 
fibre  and  matrix  displacement,  u\  =  viz  =  4)  =  itJiz  =  4)  at  the  boundaiy  of  the  process 
zone  minus  the  displacement  expected  if  the  process  zone  material  were  elastic; 

11  =  111  -  fi4 


=  -^arj  e 
k 

kE^  ^ 


(11) 


Equation  (1 1)  constitutes  the  traction  law  for  the  case  of  linearly  rising  loads. 

The  limit  of  very  fast  loading  corresponds  to  4  — >  oo,  whereupon  a  — >  0  and  17  ->  1,  since 
the  first  and  third  terms  in  Eq.  (9c)  become  negligible.  The  disturbance  then  propagates 
at  the  bar  wave  speed  in  the  matrix. 

Static  loading  is  represented  by  the  limit  4  — >  0  or  a  — >  00,  for  which  the  first  two  terms 
of  Eq.  (9c)  become  small  and  one  has  the  asymptotic  solution 


77-^ 


2afi 


For  this  limit,  substituting  Eq.  (12)  into  Eqs.  (6a)  and  (11)  yields 

/  2r 


W  11^®'^  = 


1-/ 

,  / 


172 


4E,r 


(12) 


(13a) 


(13b) 
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which  coincide  with  the  results  obtained  by  McCartney  (1987).^ 

Equations  (10)  and  (13b)  show  that  for  loads  that  increase  linearly  with  time,  the  form  of 
the  traction  law  is  identical  in  the  static  and  dynamic  cases  and 

=  •  (14) 

Thus  the  strength  of  the  effects  arising  from  inertia  is  measured  by  the  degree  to  which 
the  limit  of  Eq.  (12)  is  not  approached.  This  can  be  conveniently  summarised  by 
comparing  the  product  laPr]  to  unity. 

Large  Asymptotic  analysis  of  Eq.  (9c)  shows  that  in  the  limit  >0  la^rj  -» 

-(apf  (second,  third,  and  fourth  terms  of  Eq.  (9c)  dominant),  which 

is  a  function  of  and  not  of  ^  separately;  and  Fig.  3a  shows  that  this  limit  is  approached 
quite  closely  for  ^  >  1,  which  is  expected,  for  example,  for  composites  containing 
relatively  stiff  fibres  {Et  >  E^.  Fig.  3a  also  suggests  that,  as  an  engineering  estimate, 
inertial  effects  are  large  when 


Small  p.  The  limit  of  a  relatively  stiff  matrix  can  be  analysed  by  observing  that  a  oc  E'^ 
while  p  oc  E^ ,  and  that  Eq.  (9c)  in  the  limit  Em  oo  therefore  yields  the  limiting 

(first  and  last  terms  of  Eq.  (9c)  dominant).  Thus  lapij  -> 
(2a^  p  when  ap  ->  0.  Figure  3b  shows  this  approximation  for  three  values  of  yff  ^  1 . 
While  the  limit  is  approached  only  for  ap  so  small  that  laprj  is  also  small,  it  can 
nevertheless  be  used  as  the  basis  for  engineering  estimates  of  the  condition  for  significant 
inertial  effects,  i.e.,  the  first  significant  departure  of  laPrj  from  unity.  The  construction 
of  Fig.  3b  suggests  that  inertial  effects  will  be  significant  when  aP<  A{aP)\,  where  {ap)\ 
IS  the  value  of  the  product  ap&X  which  the  equation  {2ap)^'^p'^  =  1  is  satisfied;  i.e.. 


Marshall  et  al.  (1985)  derived  expressions  witli  different  coefficients,  but  those  of  McCartney  are  to  be 

pre  erre  since  they  correspond  to  a  proper  definition  of  the  crack  displacement  and  are  consistent  with 
conserx'ation  of  energy. 
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(17) 

This  leads  to  the  criterion  for  significant  inertial  effects  for  composites  with  relatively 
stiff  matrices  that 

{p<\)  .  (18) 

L  /  J  ^ 

The  criteria  of  Eqs.  (16)  and  (18)  coincide  when  P=\. 

Since  lapr]  <  1  always,  inertial  effects  increase  the  stiffness,  dpidii,  of  the  traction  law 
for  loads  that  increase  linearly  in  time. 

5.  Fibre  Pullout  from  a  Rigid  Matrix 

Analytical  results  have  not  been  obtained  for  other  loading  histories  for  the  composite 
problem.  (The  matrix  and  fibre  motions  evolve  in  quite  complicated  ways  for  other 
cases.)  However,  results  for  the  pullout  of  a  single  fibre  from  a  rigid  matrix  to  which  it  is 
coupled  by  friction  can  be  found  for  step  loads  as  well  as  linearly  increasing  loads.  The 
problem  considered  for  a  rigid  matrix  is  as  shown  in  Fig.  4.  The  fibre  is  loaded  on  the 
fracture  plane  by  tractions,  T(i),  which  can  be  represented  by  the  boundary  condition, 
eo{t),  for  the  axial  strain  in  the  fibre  at  z  =  0.  The  response  of  the  system  that  is  of  interest 
is  wholly  represented  by  the  load  point  displacement,  i.e.,  the  displacement,  i/o(/),  of  the 
fibre  at  z  =  0,  since  the  matrix  is  rigid  (motionless).  Specifying  a  boundary  condition  on 
the  fracture  plane,  rather  than  at  the  end  of  the  slip  zone,  as  for  the  composite  problem,  is 
preferred  here  because  it  allows  a  simpler  statement  of  the  step  loading  case. 

Full  solutions  of  this  problem  for  step  and  linearly  increasing  loads  are  given  in  Appendix 
A.  For  the  present  discussion,  the  most  interesting  features  are  the  following. 

While  the  load  point  displacement  is  reduced  by  inertial  effects  for  linearly  increasing 
loads,  for  a  step  load  it  is  increased.  For  a  linearly  increasing  load,  sa  =  kt. 
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with  the  subscript  “r”  referring  to  the  matrix  being  rigid.  The  ratio  of  Eq.  (19)  is  plotted 
in  Fig.  5.  It  is  always  less  than  unity.  In  contrast,  for  a  step  function  load  £o{t)  =  £b,  a 
constant  for  t  >  0,  all  motion  ceases  when  /  =  sqU,  at  which  point  «o  =  :  inertial 

effects  double  the  displacement  expected  from  loading  statically  to  the  same  applied  load 
Correspondingly,  the  strain  gradient  left  in  the  fibre  following  dynamic  (step)  loading  is 
exactly  half  that  found  after  static  loading. 


From  these  results,  one  can  see  that  the  following  behaviour  should  be  expected  in  the 
composite  problem  (with  a  compliant  matrix)  if  a  linearly  increasing  load  is  followed  by 
a  period  of  constant  loading.  If  the  composite  strain  has  been  increased  linearly  to  some 
value  szx  time  t\,  the  matrix  and  fibre  will  both  be  moving  at  the  boundary  of  the  process 
zone  with  (positive)  velocities  given  by  Eqs.  (10a)  and  (10b)  evaluated  at  z  =  U.  To 
achieve  this  state,  the  still-intact  composite  itself  must  be  accelerating  towards  positive  z, 
maintaining  the  same  velocity  at  z  =  4.  If  the  composite  strain  is  then  fixed  (no  longer 
increasing),  the  composite  will  slow  down  and  the  matrix  and  fibre  will  begin  to 
compress  into  the  boundary  of  the  process  zone,  under  the  influence  of  their  inertia.  This 
gives  rise  to  quite  complex  motions,  but  the  end  result  will  be  that  the  matrix  will  have 
displaced  further  relative  to  the  fibre  at  the  crack  plane  and  the  effective  crack 
displacement,  u,  will  have  increased  from  its  value  at  t  =  t\.  Analogy  with  the  problem  of 
the  fibre  being  pulled  out  of  a  rigid  matrix  suggests  that  u  might  finally  be  larger  than  if 
the  composite  strain  increased  to  the  same  level  statically.  However,  this  conjecture  must 
be  substantiated  by  numerical  solutions. 

6.  Implications  for  Dynamic  Bridged  Cracks 

Whether  or  not  dynamic  effects  in  the  pullout  process  will  be  important  in  the  problem  of 
a  propagating  bridged  crack  will  depend  on  the  loading  rate  for  the  bridging  element 
(stitch  or  rod).  The  loading  rate  will  depend  on  the  crack  profile,  which  will  depend  on 
the  specimen  shape  and  the  loading  configuration.  It  will  also  depend  on  the  crack 
propagation  rate,  but  the  onset  of  significant  dynamic  effects  can  be  evaluated  by 
considering  the  rate  of  change  of  displacements  implied  by  static  solutions  to  the 
instantaneous  load. 

The  likely  magnitude  of  inertial  effects  is  analysed  for  two  crack  problems  in  the 
o  lowmg.  To  derive  analytical  estimates  based  on  the  results  obtained  above,  the 
composite  strain  at  the  boundary  of  the  process  zone  must  be  assumed  to  rise  linearly  in 
time.  Of  course,  this  may  not  be  the  case  -  the  dynamic  crack  propagation  problem  must 
e  so  ved  with  a  self-consistently  derived  dynamic  traction  law  before  the  form  of  the 
rate  of  increase  can  be  known.  But  assuming  a  linear  increase  will  yield  insight  into  the 
likely  order  of  magnitude  of  inertial  effects. 


6.1  Steady-State,  Fully-Bridged  Crack  -  the  ACK  Limit 

Consider  first  a  crack  that  would  propagate  under  static  conditions  to  the  so-called  ACK 
imit,  m  which  a  uniform  applied  stress  comes  into  equilibrium  with  the  bridging  stresses 
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provided  by  intact  fibres  in  the  far  crack  wake  (Fig.  6).  The  ACK  limit  can  be  attained, 
for  example,  by  a  matrix  crack  in  a  ceramic  matrix  composite  in  an  infinite  specimen 
loaded  in  remote  tension  (Aveston,  Cooper,  and  Kelly,  1971;  Cox  and  Marshall,  1994). 
The  static  applied  stress,  oack,  required  to  continue  crack  propagation  in  the  ACK  limit 
is  invariant  and  the  crack  opening  displacement  in  the  far  crack  wake  is  uniform  (Fig.  6). 
For  a  composite  in  which  the  bridging  elements  are  coupled  to  the  matrix  by  uniform 
friction,  static  analysis  gives  the  following  well-known  results  (e.g.,  McCartney,  1987; 
Appendix  B  in  Massabo  and  Cox,  1999;  Marshall  and  Cox,  1988).  The  critical  stress, 
Oack,  is  given  by 


^ACK  - 


(21) 


where 


r 


p  _  IfE  Eft 


(22) 


and  Gc  is  the  critical  crack  tip  energy  release  rate  for  the  matrix  crack.  The  bridging 
stress  approaches  close  to  oxck  at  a  characteristic  distance,  /ack,  from  the  crack  tip  given 
by 


'ACK 


kE 

T 


(23) 


where  E  is  an  elastic  constant  that  depends  on  the  degree  of  anisotropy  of  the  composite 
(defined,  e.g.,  in  Cox  and  Marshall,  1991).  If  the  composite  strain  on  the  boundary  of 
the  process  zone  in  the  wake  of  the  ACK  crack  is  assumed  to  rise  linearly  in  time,  as  in 
Eq.  (8),  then 


(24) 

^ACK  /^d 

where  Cd  is  the  velocity  of  propagation  of  the  delamination  crack  tip,  which  is  assumed 
constant.  Using  Eq.  (9d),  the  product  afS,  which  determines  the  magnitude  of  inertial 
effects  for  a  linearly  rising  load,  is  then 

(ACK  crack)  .  (25) 

\6  f  E  Ef  Cj 
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For  the  square  root  traction  law,  p  oc  tt^'\  the  friction  stress,  r,  the  fibre  radius,  R,  and  the 
critical  energy  release  rate  for  the  matrix,  Gc,  do  not  appear  in  this  expression,  although 
they  appear  in  /ack  and  oxck  separately  (Eqs.  (21)  and  (23)). 

For  composites  satisfying  p  >  \,  the  criterion  for  significant  inertial  effects  in  the 
bridging  mechanism  is  that  ap<2  (Eq.  (15)).  For  elastically  homogeneous  composites, 
E^E  and  aP  depends  only  on /and  the  ratio  of  the  wave  speed  and  the  crack 
velocity  a.  For  example,  for/=  0.5,  will  have  a  value  less  than  2  if  cj  >  (;z/32)Cm. 
The  velocity  of  a  crack  propagating  in  the  ACK  configuration  has  not  previously  been 
calculated.  However,  if  the  crack  is  fully  bridged,  as  in  an  unnotched  composite 
exhibiting  multiple  matrix  cracking,  crack  growth  is  known  to  be  unstable  (e.g.,  Cox  and 

Marshall,  1994),  so  that  crack  velocities  satisfying  cj  >  (/il32)cm  would  appear  to  be 
feasible. 

For  composites  satisfying  /?  <  1  (/  ~  0. 1  would  appear  to  be  easily  attainable),  inertial 
effects  will  appear  at  crack  velocities  lower  by  a  factor  of (Eq.  (18)). 

Thus  inertial  effects  are  predicted  to  be  significant  for  the  ACK  crack  configuration  for  at 
least  some  common  composites. 

6.2  W edge-Loaded  Double  Cantilever  Beam 

The  double  cantilever  beam  (DCB)  specimen  loaded  dynamically  by  a  flying  wedge 
offers  a  relatively  simple  experimental  approach  to  the  mode  I  dynamic  delamination 
problem  (Fig.  7).  The  test  is  especially  attractive  for  studying  the  bridging  effects 
supplied  by  through-thickness  reinforcement  (e.g.,  stitches  or  rods)  in  laminates.  An 
estimate  of  the  likely  role  of  inertia  in  bridging  by  stitches  or  rods  follows.  In  the  context 

of  this  paper,  the  “fibre”  refers  in  this  case  to  a  stitch  or  rod  and  the  “matrix”  to  the 
laminate. 

For  an  increasing  bridging  traction  law,  i.e.,  a  law  for  which  dp/dit  >  0,  the  crack  surfaces 
at  the  crack  tip  are  predicted  to  come  into  contact  in  the  DCB  specimen  under  static 
loading  at  a  certain  crack  length  (Massabo  and  Cox,  2000).  The  crack  will  be  arrested  at 
this  point,  since  the^  crack  tip  energy  release  rate  must  then  vanish.  Furthermore,  the 
crack  surfaces  remain  in  contact  when  the  applied  load  is  increased  so  that  the  crack 
remains  arrested  and  failure  eventually  ensues  by  another  mechanism.  An  analogue  of 
t  IS  interesting  phenomenon  of  crack  tip  closure  and  crack  arrest,  which  is  a  result  of 
bridging  effects,  might  also  be  anticipated  in  dynamic  loading  (Beyerlein  et 
a ).  An  interesting  case  therefore  is  to  consider  the  rate  of  opening  of  a  crack  that 
h&s  Errcstcd  snd  is  being  loaded  further  by  the  flying  wedge. 

Simple  analytical  expressions  for  the  opening  displacements  are  available  in  the  static 
case  for  a  linear  bridging  law, /?  =  Pju,  and  a  composite  with  moderate  levels  of  elastic 
anisotropy.  If  /  is  suitably  large  in  Fig.  7,  then  the  bending  moment  at  the  notch  root,  .v  = 
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0,  will  dominate  over  the  shear  load  at  x  =  0.  In  this  case,  the  crack  will  arrest  after 
propagating  a  distance  a\  given  by  (Massabo  and  Cox,  2000) 

(26) 


where 


A 


(27) 


with  -  Vxy  Vyx)  (plane  strain  conditions  assumed)  and  h  =  12//?^,  with  h  the 

laminate  half-thickness  (see  Fig.  7  for  coordinates  and  dimensions).  Insight  into  the 
magnitude  of  dynamic  effects  in  the  bridging  phenomenon  can  be  gained  by  considering 
the  rate  at  which  the  crack  opens  at  the  end  of  the  bridging  zone  (x  =  0  in  Fig.  7)  under 
the  influence  of  the  wedge  load  when  the  crack  length  is  fixed  aX  a  =  a\.  From  the  results 
of  Massabo  and  Cox  (2000),  the  load  point  displacement,  5,  at  the  point  where  the  wedge 
contacts  the  specimen  can  be  related  for  large  /  to  the  crack  opening  displacement,  Uq  ,  at 
X  =  0  (the  last  intact  bridging  element)  by  (see  Appendix  B)^ 


S  =  -^ 


Ur, 


(28) 


The  composite  strain  at  the  boundary  of  the  process  zone  at  x  =  0  is 


_ P^Uq 

C  — 


(29) 


while,  for  a  wedge  subtending  an  angle  2^  and  moving  at  velocity  v, 


w, 


dS 

—  =  v^X?in<j) 


(30) 


If  Vw  is  constant,  Eqs.  (28  -  30)  yield  a  linear  rate  of  increase  for  s,  where,  in  the  notation 
of  Eq.  (8), 


^  Tliese  estimates  are  based  on  results  for  a  linear  bridging  law.  Equivalent  results  for  a  quadratic  law,  p  oc 
which  is  expected  for  static  loading  of  fibers  coupled  to  the  matrix  by  friction,  have  not  been 
published.  However,  the  estimated  times  for  reaching  a  given  displacement  and  composite  strain  are  likely 
to  be  of  the  same  order  of  magnitude  as  for  a  linear  law. 
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k 


A 

1 — 
-s: 

2 

- 1 

_ 1 

tan^ 

~~F~ 


and  thence 


(31) 


ap  = 


— _ ^  ^  -^3 

V3  1-/  ^E'^P^h  hR  tan^zJ 


(wedge-loaded  DCB)  . 


(32) 


Here  refers  to  the  bar  wave  speed  in  the  laminate  (“matrix”)  in  the  through-thickness 
direction.  The  modulus^  is  that  of  the  composite  (with  stitches  or  rods)  in  the  through- 
thickness  direction,  while  Em  is  that  of  the  laminate  (without  stitches  or  rods)  in  the  same 
direction.  If/is  small,  which  is  the  usual  case,  Em  «  E^,  and  the  factor  //(l  -/)  with 

/=  0.05  a  typical  value.  Experimental  data  for  typical  stitched  laminates  show  «  100 
MPa/mm  (Massabo  et  al.,  1998;  Turrettini,  1996).  Taking  E'^  =  60  GPa  and  ^  =  6  mm, 

the  factor  has  the  representative  value  6  GPa.  The  ratio  t/^E'^p^h  will  have 

typical  values  ~  lO'^  -  lO'^  (r=  6  -  60  MPa).  With  /  =  20  mm,  i?  =  1  mm,  and  <!>  =  20°, 
for  example,  Eq.  (32)  reduces  to  the  particular  numerical  estimate 


0.01^  <ay9<  0.1^ 
v,„  V... 


For  composites  for  which  p  >  1  (including  most  stitched  laminates),  inertial  effects  will 
be  significant  (i.e.,  ap  <  2)  provided  Vw  >  cJ200  (r  =  6  MPa)  or  Vw  >  cJ20  (t  =  60 


Such  wedge  velocities  are  eminently  attainable  and  inertial  effects  in  the  bridging  law  are 
likely  to  be  significant  in  the  wedge-loaded  DCB  test. 

7.  Discussion 

The  solutions  presented  here  are  based  on  the  approximation  that  displacements  are 
unctions  of  z  only.  For  static  loading,  this  is  an  accurate  approximation  as  long  as  the 
slip  distance  is  large  compared  to  the  fibre  diameter  (Hutchinson  and  Jensen,  1990).  The 
condition  that  the  slip  zone  length  should  be  much  larger  than  the  fibre  diameter  will  be 

satisfied  if  the  friction  stress,  r,  takes  values  typical  in  brittle  matrix  composites  or 
polymer  textile  composites. 

For  dynamic  problems,  the  conditions  under  which  solutions  that  depend  on  z  only  will 
e  accurate  have  not  yet  been  studied.  The  solutions  presented  in  this  work  will  at  least 

provi  e  tentative  limits  against  which  three-dimensional  numerical  solutions  can  be 
compared. 
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Whether  the  approximation  that  displacements  depend  on  z  alone  is  accurate  will  also 
depend  on  the  presence  of  a  nonzero  debond  energy  (energy  required  for  separation  of  the 
fibres  and  the  matrix).  Here  the  debond  energy  was  assumed  to  be  zero  (fibres  and 
matrix  already  chemically  separated).  Dynamic  analysis  of  materials  in  which  the 
debond  energy  is  not  small  will  be  significantly  more  complicated. 

For  stitched  laminates,  one  of  the  cases  assessed  in  Section  5,  the  requirement  that  the 
slip  or  process  zone,  4,  be  large  compared  to  the  fibre  (stitch)  diameter  (typically  1  mm) 
implies  that  it  will  also  be  comparable  to  or  greater  than  a  typical  laminate  half-thickness, 
/?  (~  6  mm).  When  the  value  of  4  predicted  for  an  infinite  body  would  exceed  h,  stress 
waves  will  reflect  from  the  laminate  surface,  complicating  the  micromechanics  of  the 
process  zone.  These  effects  are  not  studied  here. 

A  feature  of  the  wave  solutions  for  loads  that  increase  monotonically  from  zero  with  time 
is  that  stress  disturbances  propagate  away  from  the  fracture  plane  at  velocities  other  than 
either  the  matrix  or  fibre  bar  wave  speed.  This  is  the  case  for  both  the  composite  problem 
(Section  4)  and  the  problem  of  a  fibre  being  pulled  out  of  a  rigid  matrix  (Nikitin  and 
Tyurekhodgaev,  1990;  Appendix  A).  In  the  composite  problem,  the  front  velocity  is 
always  less  than  the  matrix  bar  wave  speed.  The  possibility  exists  that  the  front  velocity 
will  exceed  the  bar  wave  speed  in  the  fibres,  but  numerical  checks  show  that  this  can 
happen  only  for  very  high  (and  unlikely)  ratios  of  fibre  to  matrix  density. 

No  account  has  been  attempted  here  of  the  conditions  under  which  the  strain  at  the 
boundary  of  the  process  zone,  s,  can  be  expected  to  rise  linearly  in  time,  apart  from  the 
special  case  of  a  wedge  loaded  DCB  specimen.  The  relationship  between  e  and  the  far 
field  conditions  in  a  dynamically  loaded  body  will  generally  be  complicated  and  revealed 
only  by  computational  solutions  of  the  dynamic  stress  propagation  problem  for  the  whole 
body,  including  the  process  zone.  The  particular  problem  of  finding  far  field  loading 
conditions  that  will  result  in  £(/)  being  linear  in  time  is  very  challenging. 

8.  Conclusions 

Some  analytical  results  have  been  presented  for  the  problem  of  bridging  by  the 
mechanism  of  fibre  pullout  when  the  inertia  of  the  fibre  and  the  matrix  are  taken  into 
account.  Simple  criteria  have  been  specified  for  significant  inertial  effects  in  the  bridging 
mechanism  in  representative  mode  I  crack  propagation  problems.  Inertial  effects  in  the 
bridging  mechanism  will  often  be  significant  for  a  matrix  crack  propagating  dynamically 
in  the  steady  state  ACK  limit  in  a  brittle  matrix  continuous  fibre  composite  or 
delamination  in  laminates  reinforced  by  through-thickness  stitching  or  rods. 

Significant  inertial  effects  are  favoured  by  low  fibre  volume  fraction,  low  friction  stress, 
low  matrix  bar  wave  speed,  and  low  fibre  modulus;  and  high  fibre  diameter  and  high 
matrix  modulus.  If  the  fibre  modulus  is  high  enough  (relative  to  the  matrix  modulus),  the 
criterion  for  inertial  effects  becomes  independent  of  fibre  modulus.  The  matrix  density 
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enters  the  criterion  only  through  the  matrix  bar  wave  speed.  The  criterion  is  always 
independent  of  the  fibre  density. 

For  pullout  or  bridging  stresses  that  rise  linearly  in  time,  the  instantaneous  crack 
displacement  is  less  in  the  presence  of  inertial  effects  than  it  would  be  under  static 
loading  to  the  same  bridging  stress.  However,  solutions  for  pullout  from  a  rigid  matrix 
suggest  that,  if  the  bridging  stress  rises  rapidly  and  is  then  held  at  a  constant  value,  the 
crack  displacement  when  all  particle  motion  finally  stops  will  be  greater  than  it  would 
have  been  under  static  loading  to  that  stress  level.  Thus  regimes  of  both  hardening  and 

softening  of  the  bridging  traction  law  due  to  inertial  effects  can  be  expected  in  bridged 
crack  problems. 
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Appendix  A.  Fibre  Pullout  from  a  Rigid  Matrix. 

Solutions  are  exhibited  here  to  the  problem  described  in  Fig.  4  and  related  text.'^  Since 
the  matrix  is  rigid,  it  experiences  no  motion.  The  motion  of  the  fibre  is  governed  by  the 
wave  equation  of  Eq.  (2a).  The  problem  is  defined  by  stating  the  value  of  the  load 
applied  to  the  end  of  the  fibre  at  z  =  0,  rather  than  the  matrix  strain  at  the  end  of  the 
process  zone. 

Boundary  Conditions  for  Front  of  Propagating  Disturbance 

Let  the  location  of  the  front  limiting  the  extent  of  the  stress  disturbance  caused  by  waves 
propagating  along  the  fibre  be  written  U  =  tic^,  where  the  applied  load  is  turned  on  at 
time  /  =  0.  (The  subscript  “r”  is  used  for  this  problem  as  a  mnemonic  for  “rigid”  matrix.) 
The  jump  in  the  stress,  A  of,  and  jump  in  particle  velocity,  Ai/f  (where  =du^ldt), 

across  the  front,  i.e.,  from  z  >  rjc^  to  z  <  Tjcf  must  satisfy  the  energy-conserving  relation 
expected  from  integration  of  the  impulse  across  the  front,  namely  (see,  e.g.,  Achenbach, 
1973) 


Ai/f 


Acf 

Pf^O 


(A.1) 


where  vq  is  the  velocity  of  the  front,  vq  =  C[8{T}t)l8t .  To  ensure  integrity  of  the  fibre,  the 
displacement  must  be  continuous  across  the  front,  i.e., 

«f(z,r)  =  0  .  (A.2) 

Two  front  conditions  may  now  be  distinguished.  If  a  stress  discontinuity  exists  at  the 
front,  Aaf  >0,  then  kinematic  considerations  along  with  Eq.  (A.1)  necessitate  that  (see, 
e.g.,  Achenbach,  1973) 

(Acrf>0)  (A.3a) 


or,  equivalently,  vq  =  Cf.  If  the  stress  is  continuous  across  the  front,  i.e.,  af  =  0  at  z  =  Zf, 
then  the  velocity  of  the  front  remains  indeterminate; 


0<  7<  1 


(AcTf  =0) 


(A.3b) 


4 

Tlie  same  problem  of  an  end-loaded  rod  damped  by  friction  was  also  solved  for  various  cases  by  Nikitin 
and  colleagues  in  the  1960’s.  These  slightly  arcane  Russian  publications  are  summarised  in  Nikitin  and 
T>mrekhodgaev  (1990),  which  itself  has  been  cited  only  once  previously  in  the  western  literature.  In  the 
following,  dynamic  solutions  are  developed  for  tlie  cases  of  particular  interest  to  the  delamination  problem, 
with  emphasis  on  tlie  relationship  betw^een  tlie  dynamic  results  and  the  static  limit.  Some  minor  errors  in 
Nikitin  and  Tjiirekhodgaev  for  the  case  of  loading  that  is  linear  in  time  are  corrected. 
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Solutions  obeying  both  Eqs.  (A.3a)  and  (A.3b)  will  be  demonstrated  in  the  following. 
Step  function  Load 

Consider  first  a  step  function  load  such  that 


-EfS^ 


for  /  >  0 


(A.4) 


Since  a  stress  discontinuity  is  implied  by  a  propagating  step  load,  the  front  conditions 
given  by  Eqs.  (A.l),  (A.2),  and  (A.3a)  are  expected  to  apply.  The  solution  to  the  wave 
equation,  Eq.  (2a),  that  satisfies  these  boundary  conditions  is 


with  Tj  1  (front  velocity  vq  —  Cf)  always.  The  particle  velocity  is  thus 


RE, 


t  £qC^  —  Cj* 


(A.  6a) 


where 


t. 


EfR 

TC[ 


(A.6b) 


and  is  independent  of  z:  within  the  sliding  zone  0  <  z  <  the  fibre  moves  with  uniform 
velocity.  The  particle  acceleration  is  given  by 


5^  r_  2  _£f 
dt^  REf  "  t^ 


(A.7) 


which  is  uniform  in  time.  Motion  stops  when 


t=£btr 


(A.8) 


This  is  the  characteristic  time  of  the  dynamic  process  for  the  case  of  a  rigid  matrix.  At  /  = 
cotr,  the  sliding  zone  has  advanced  a  distance 


;  _  .  EcR 

T 


{t=SoU)  (A.9) 
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The  Strain  along  the  rod  has  the  distribution 

ditc  T 

dz  “  EEi 

The  load  point  displacement,  r/o,  (i.e.,  the  displacement 
according  to 


(A.10) 

at  z  =  0)  evolves  in  time 


Wq 


l- 


1  r 

2so 


(A.  11) 


When  t  =  soti. 


E  R 

"o  =  ^^0^  (/  =  SoU)  (A.  12) 

which  is  exactly  twice  the  displacement,  developed  under  static  loading  (see 

below).  The  final  slip  length,  4,  is  also  twice  the  slip  length  for  static  loading;  and  the 
strain  gradient,  du^/dz ,  is  half  that  developed  under  static  loading. 

One  important  consequence  of  the  last  fact  is  that  when  motion  is  arrested  at  t  =  no 
further  stress  relaxation  is  required:  the  stress  gradient  at  t  =  sot,  is  only  half  the 
maximum  gradient  that  can  be  supported  in  the  fibre  by  the  friction  tractions.  According 
to  the  constitutive  behaviour  assumed  for  the  friction  phenomenon,  Eq.  (Ic),  the  fibre  will 
remain  motionless.  Therefore,  the  configuration  predicted  for  t  =  sqU  is  the  final 
configuration.  Since  Eq.  (A.  10)  predicts  that  o}(/r)  =  0  at  /  =  Sotj,  the  front  condition  will 
switch  at  this  instant  to  Eq.  (A.3b):  the  condition  of  Eq.  (A.  1)  leads  to  an  indeterminate 
front  velocity. 

When  the  fibre  arrests  at  /  =  soU,  the  particle  acceleration  in  the  slip  zone  falls 
instantaneously  to  zero  from  the  constant  value  of  Eq.  (A.  7).  At  the  same  instant,  the 
friction  tractions  along  the  slip  zone  fall  from  the  uniform  value,  r,  which  they  had  during 
the  motion  of  the  fibre,  to  r/2,  the  value  required  to  sustain  the  stress  gradient  in  the  fibre 
at  the  time  of  arrest.  This  discontinuous  change  in  the  friction  stress  is  a  direct 
consequence  of  the  constitutive  behavior  embodied  in  Eq.  (1). 

For  the  step  load  case,  Eq.  (A.  5)  yields  the  following  energy  analysis.  The  work  done  by 
the  load  at  time  t  <  eoU  is  given  by 


W,  =  TtR^EcSf^Ur,  = 

r '  1 

'i-i 

2 

'  11 

r 

•1 

Uo^j. 

The  kinetic  energy  in  the  moving  fibre  at  time  t  is  given  by 
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(A.13) 
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fT,  =  dz  =  _L 

2  {dt  J  2t 


(A.  14) 


The  strain  energy  in  the  fibre  at  time  t  is  given  by 


O  i_f -Ll+if _lT 

*  2  [  dz  J  2t  I  sJ.  J  31  £J. 


(A.15) 


These  energy  dissipated  in  friction  up  to  time  t  is  given  by 


nE^slE?  i  t 


JVf  =  -/gf'  iTtRnidz  = 


1-^U 


(A.  16) 


The  energy  terms  are  plotted  in  Fig.  A.l.  The  total  work  done  at  time  t  =  sqU  is  exactly 
three  times  that  done  in  static  loading  to  (see  below). 

Linearly  Increasing  Load 

Consider  next  the  case  of  an  applied  strain  that  is  given  by  the  linear  law 


>^f(0)=£f^  =E,kt  (/>0) 

&  z=0 


(A.17) 


In  this  case,  a  physically  meaningful  solution  is  found  only  for  the  front  condition  given 
by  Aof  =  0  and  Eq.  (A.3b).  With  these  boundary  conditions,  Eq.  (2a)  has  the  solution 


t/f(z,r)  =  - 


/l  +  ikt, 


T  1 

REf  1-1]' 


-{z-r]cdY 


(Z  <  TJCft) 


(A.  18) 


Tl  =  kLj  ^l\  +  {kt,Y  +1 


(A.19) 


where  U  is  the  characteristic  time  of  Eq.  (A.8).  Thus  the  particle  velocity  is  given  by 
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The  Static  limit  in  this  case  can  be  found  by  taking  the  limit  k^O.  One  finds 


2t  lx  '' 

(k  — ^  0) 

(A.24a) 

dt  \_  T  _ 

(k-^O) 

(A.24b) 

(k^O) 

(A.24c) 

4r 

(k^O) 

(A.24d) 

where  the  superscript  “(st)”  indicates  the  result  obtained  by  static  analysis  for  loading  to 
the  instantaneous  value,  ch(0,t),  of  the  applied  tractions  (see  below).  In  the  static  limit, 
Eq.  (A.24b)  also  shows  that  tj  ->0  (consistently  with  Eq.  (A.24a));  while  in  the  limit  of 
very  rapid  loading  (k  ->  oo),  or  vanishing  friction  (r  0),  Eq.  (A.24b)  leads  to  I, 
i.e.,  the  front  reverts  to  propagating  at  the  bar  wave  velocity. 
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Equations  (A.23)  and  (A.24d)  allow  the  dynamic  load  point  displacement  to  be  related 
very  simply  to  the  displacement  expected  for  static  loading  to  the  same  instantaneous 
applied  load: 


1 

o 

1 

+ 

•-I 

1 

1 

(A.25) 


The  same  characteristic  time,  t^,  identifies  cases  where  dynamic  effects  are  important  in 
both  the  step  and  linear  loading  cases. 

The  Static  Problem 


Under  static  loading  to  stress  oo  -  EiSo,  force  equilibrium  leads  to  a  linear  stress  gradient 
along  the  fibre 


^(st)  _  ^ 

Of-  -  o^ 


1- 


Ixz 

cTqR 


so  that  the  slip  length,  where  =  0,  is  given  by 


(A.26) 


/(St) 

'r  ~  «  *0 

2r 


(A.27) 


The  load  point  displacement  can  be  found  by  integrating  the  strain  along  the  fibre  and  is 
given  by 


,(st) 


"o  ^  =  - 


EfR 


At 


The  total  work  done  by  the  load  in  the  static 


=  f'rrR^adu  =  — 

6r 


problem  is  given  by 

3 


where  here  u  and  o denote  the  displacement  and  stress  at  z  =  0  during  loading. 


(A.28) 


(A.29) 
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Appendix  B.  The  Wedge-Loaded  DCB  Specimen  with  Large  Scale  Bridging. 

In  a  specimen  in  which  shear  deformation  may  be  neglected  (material  anisotropy  not  too 
large)  and  in  which  a  linear  bridging  law,  p  acts  on  the  bridged  part  of  the  crack 
(Fig.  7),  analytical  expressions  can  be  found  for  the  characteristics  of  crack  propagation 
(Massabo  and  Cox,  2000).  The  crack  displacement  profile  is  given  by 

i/  =  e*^(ciCOS&Ji:  +  C2sin6x)+e“*^(c3CosZ»x+C4sinZ>x)  (B.l) 


where 


b  = 


(B.2) 


and  Id  -  A^/12,  =  EJ(\  -  VyxKxy),  and  E^  and  Vyx  and  Vxy  are  Young’s  modulus  and 

Poisson’s  ratios  for  the  laminate,  respectively.  The  constants  Ci  are  determined  by 
boundary  conditions.  At  the  crack  tip,  ?/  =  0  and  the  bending  rotation  (f)  =  -  diildx  =  0 
(shear  deformation  neglected);  while  at  the  notch  root,  x  =  0,  one  can  write  generally  that 

//  =  Uq  and  Observing  that  the  bending  moment,  M  =M(x  =  0),  and  shear  force, 

O  =  0(x  =  0),  at  X  =  0  are  related  by 


K--EL 

Q~  P 


(B.3) 


where  P  is  the  load  applied  by  the  wedge  and,  for  a  beam  with  negligible  shear 
deformation,  0  =  -(£'x'fd)d^i//dx^  and  M  =  -(£x'/d)d^?//dx^,  one  finds  by  solving  for  the 
coefficients  ci  in  (B.l)  that 


j  ~  6  ^  + 1  +  llbip  —  l) 

- e-'-lH-ftp+l) 


(B.4) 


and  thence 


e-’^  -\+lb[e-’‘  +  \) 


(B.5) 


If  the  loading  arm  length,  /,  is  not  short,  then  the  deflection  at  the  point  of  contact  with 
the  flying  wedge,  5,  will  be  much  greater  than  hq  and  therefore  5  can  be  estimated  as  the 
deflection  expected  for  a  cantilever  beam  with  built-in  end  condition  at  x  =  0; 
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and  thus 


In  the  limit  that  lb  is  large, 


S 


l.\Ml 


(B.6) 


(B.7) 


(B.8) 
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Figure  Captions 

1.  (a)  Schematic  of  a  crack  bridged  by  fibres,  showing  process  zone  where  relative 

disp  acement  exists  between  fibres  and  matrix,  (b)  Idealisation  of  the  bridged 
crack  problem  with  process  zone  replaced  by  surface  tractions  acting  on  Ihe 
fracture  surfaces.  ^ 

2.  Schematic  of  the  dynamic  pullout  problem  in  a  representative  volume  near  the 
fracture  plane. 

3.  The  product  2 apt],  which  indicates  the  relative  importance  of  inertial  effects  in  the 
bridging  phenomenon,  (a)  Numerical  results  for  representative  values  of  B.  fb) 

umencal  results  for  la^r]  compared  with  asymptotic  limit  for  large  (dashed 

curves).  The  value  of  is  marked  for  0.01  and  is  to  be  compared  with  the 

value  of  ap  at  which  laprj  «  0.9,  which  is  taken  as  a  representative  cutoff  for 
significant  inertial  effects. 

4.  Fibre  pullout  from  a  rigid  matrix. 

5.  The  ratio  of  dynamic  and  static  load  point  displacements  for  a  fibre  pulled  out  of  a 
rigid  matrix  by  a  load  that  rises  linearly  with  time. 

6.  A  matrix  crack  propagating  in  the  ACK  limit. 

7.  A  double  cantilever  beam  specimen  loaded  by  a  flying  wedge. 

A.  1  The  distribution  of  energy  as  a  function  of  time  for  a  fibre  pulled  out  of  a  rigid 
matrix  by  a  step  load.  ° 


INERTIAL  EFFECTS  IN  PULLOUT 


/ 


H 


n/ 

: 

1  r 

ii 

1 

I  I 

j 

i  I 

■j 

■4 

- 

- 

] 

/ 

iv 

- 

i 

i. 

u 

^  crack  tip 


Boundary  of 
process  zone 


Figure  1 


7‘ 

matrix  crack 
plane 


of  matrix 


RSCCSMCmO«5C<- 


Figure  2 


INERTIAL  EFFECTS  IN  PULLOUT 


a> 


front  of  stress 
y' propagation 


frictional  tractions,  r 


n  n  0 
T(/)  = 


n 


RSC0M3Cc»  080500 


Figure  4. 


Figure  5. 


P 


71 


INERTIAL  EFFECTS  IN  PULLOUT 


bridged 

zone 


RSCOS?!  Co*  072500 


Figure  7. 


energy 


rca^R^ 

-tE 


time,  //fj 


RSCC574COZC80500 


Figure  A.l 


72 


Simple  form  for  the  Dynamic  G-K 
relationship  in  Orthotropic  Materials 

N.  Sridhar  and  B.  N.  Cox 

Rockwell  Science  Center 
1049  Camino  Dos  Rios 
Thousand  Oaks,  CA  91360 
U.S.A. 

and 

1.  J.  Beyerlein 

Los  Alamos  National  Laboratory 
Los  Alamos,  New  Mexico 
U.SA 


To  be  submitted  as  a  “brief  note”  to  J.  Applied  Mechanics 
_ February,  2001 


73 


1.  Introduction 


This  note  presents  explicit  relations  between  the  dynamic  crack  energy  release  rate 
and  crack  tip  stress  intensity  factors  in  a  form  most  amenable  to  computational  work  for  a 
crack  propagating  along  any  principal  axis  in  an  orthotropic  material.  The  results 
presented  in  this  form  may  be  particularly  useful  for  work  on  delamination  cracks  in 
laminated  composites,  which  are  orthotropic  in  many  important  applications. 

Prior  work  in  the  literature  that  treats  the  problem  of  interest  here  include  the  analysis 
of  a  mode  I  crack  in  an  orthotropic  medium  by  Piva  and  Viola  (1988);  and  the  analysis  of 
a  mixed  mode  crack  in  a  transversely  isotropic  medium  by  Wu  (1989).  Yang  and  co¬ 
workers  (1991)  presented  results  for  a  mixed  mode  crack  in  an  orthotropic  material,  but 
in  a  form  that  is  not  best  suited  for  immediate  computation.  By  an  adaptation  of  the  Stroh 
method  of  analysis  used  by  Wu,  the  results  of  Yang  et  al.  are  rendered  here  in  a  different 
form  better  suited  to  numerical  evaluation.  The  role  of  orthotropic  anisotropy  in 
determining  the  velocity  dependence  of  Gd  is  then  illustrated  for  delamination  cracks  in 
typical  ceramic,  metal  and  polymer  matrix  laminated  composites. 


2.  Results 

The  relation  between  the  energy  release  rate  and  the  crack  tip  stress  intensity  factor  is 
central  to  fracture  analysis.  It  appears,  for  example,  in  the  derivation  of  stress  intensity 
factors  from  applications  of  the  J-integral  and  in  the  derivation  via  energy  arguments  of 
integral  equations  for  bridged  cracks  from  weight  functions.  Details  of  the  latter  for 
dynamic  cracks  with  large  scale  bridging  in  the  crack  wake  will  appear  elsewhere. 

Freund  (1972)  showed  that  the  dynamic  crack  energy  release  rate  for  an  extending 
crack  in  an  elastic  body  can  be  written  as  the  modified  J-integral,  which  leads  to  the  G-k 
relationship  in  the  following,  convenient  form: 

Gd  =  (1) 
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where  Gd  is  the  dynamic  crack  energy  release  rate  and  k(t)  is  the  instantaneous  crack  tip 
stress  intensity  factors  (see  Eq.  4.19  in  Wu  and  references  therein).  The  elements  of  the 
L'*  matrix  are  universal  functions,  in  the  sense  that  they  are  independent  of  the  details  of 
the  applied  loading  or  on  the  configuration  of  the  body  being  analyzed.  We  present 
results  for  the  L  matrix  when  the  material  possesses  orthotropic  symmetry  and  when  the 
crack  is  extending  along  one  of  the  principal  axis.  In  this  case,  the  off-diagonal  elements 
of  the  matrix  are  zero  and  the  diagonal  elements  of  the  L’*  matrix  is  just  dependent  on  the 
elastic  properties  of  the  orthotropic  medium  and  the  instantaneous  value  of  the  crack  tip 
velocity.  For  an  orthotropic  material  with  the  crack  propagating  along  the  1  direction  (and 
direction  3  being  the  plane  strain  direction),  elements  of  the  matrix  are; 
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where  the  y’s  are 
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(3) 

where  Q  are  coefficients  of  the  stiffness  matrix,  p  is  the  density  of  the  medium  and  v  is 
the  instantaneous  velocity  of  the  crack  tip.  The  expressions  for  P  are: 
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= 

(4) 

where 

r,  =C,,C22  -C,^2  -(C22  +Q6)/7v^ 

and 

/'2  =  4C22Q  (C, ,  -  pv^  )(Q  -  )  . 

The  L'S,3  matrix  elemenet  is 

(5) 

We  can  also  define  the  generalized  Rayleigh  wave  function  R(v)  as 

Riy)  =  rl-nVi  (6) 

and  the  Rayleigh  wave  speed  (vr)  is  obtained  by  setting  R(vr)  =  0. 

In  the  next  section,  the  diagonal  elements  L''i,i  and  u'2.2  are  plotted  numerically  for 
cracks  propagating  along  different  principal  axes  in  various  representative  composites. 
The  L*'  matrix  elements  are  functions  of  the  elastic  properties  and  the  instantaneous  crack 
velocity.  Each  diagonal  element  approaches  the  corresponding  static  value  as  v  — ^  0  and 
has  the  property  of  0[(v^  —  v)“*  ]  as  v  ->  ,  where  Vr  is  the  Rayleigh  wave  speed. 

3.  Discussion 

In  this  section,  we  will  present  results  for  the  the  mode  I  and  mode  II  contributions  to 
the  crack  energy  release  rate,  as  represented  by  the  and  L‘^2.2  matrix  elements,  for 
various  representative  composites  that  possess  orthotropic  symmetry  and  where  the  crack 
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is  extending  along  one  of  the  principal  axis.  These  results  are  presented  graphically  as  a 
function  of  the  orthotropy  parameters  and  the  instantaneous  crack  tip  velocity.  The  mode 
III  contribution  to  the  crack  energy  release  rate  is  a  simple  expression  (see  Eq.  5)  and 
hence  will  not  be  examined  further. 

In  the  results  presented  below,  the  orientation  is  such  that  the  principal  axes  of 
orthotropic  symmetry  are  aligned  with  the  coordinate  axis,  the  crack  lies  in  the  1-3  plane 
and  the  crack  is  propagating  along  the  1  direction  (with  direction  3  being  the  plane  strain 
direction).  Table  I  shows  the  stiffness  tensor  for  various  representative  composites,  the 
corresponding  orthotropy  paramaters  and  the  Rayleigh  velocity  (vr)  normalized  by  the 
shear  wave  velocity  (vs).  Also  examined  in  Table  I  are  cases  where  the  crack  is  running 
along  the  different  principal  directions  of  orthotropic  symmetry.  The  orthotropy 
parameters  X  and  p  as  presented  by  Suo,  et  al.  (1991)  are: 


and 


P  = 


*12  +  2*^ 


where 


S^j  {plane  stress) 

~  '  C  t  ■  \ 

Oy - r — {plane  strain) 

*33 

and  where  Sy  are  the  elements  of  the  compliance  tensor. 

Figure  1  shows  the  variation  of  the  normalized  L'',.,  /  ,  where  [L-',.,]s,a,ic 

is  the  value  of  L'u  when  v— >  0,  as  a  function  of  the  normalized  crack  velocity  v/vr , 
where  Vr  is  the  Rayleigh  wave  speed.  This  variation  is  shown  for  the  different  cases 
listed  in  Table  I.  However,  it  is  clear  from  the  figure  that  the  variation  in  the  crack 
energy  release  rate  can  be  rationalized  in  terms  of  the  orthotropy  parameters  X  and  p.  For 
the  same  normalized  crack  speed,  the  normalized  mode  I  contribution  (L'’,.,  /  [L'S.iJstatic)  to 
the  dynamic  crack  energy  release  rate  increases  as  X  monotonically  increases  and  as 
p  monotonically  decreases. 

Figure  2  similarly  shows  the  variation  of  the  normalized  v'2,2  /  [L‘'2,2]static ,  where 
[L  2.2]statio  is  the  value  of  L  '2,2  when  v  — >  0 ,  as  a  function  of  the  normalized  crack  velocity 

v/vr,  where  Vr  is  the  Rayleigh  wave  speed.  This  variation  is  shown  for  the  different  cases 
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listed  in  Table  I.  As  in  the  mode  I  case,  the  variation  in  the  crack  energy  release  rate  can 
be  rationalized  in  terms  of  the  orthotropy  parameters  X  and  p.  For  the  same  normalized 
crack  speed,  the  normalized  mode  II  contribution  (L’'^/  [L'^ijstatic)  to  the  dynamic  crack 
energy  release  rate  increases  as  X  increases  and  as  p  decreases.  However,  as  opposed  to 
the  mode  I  case,  the  variation  on  X  and  p  is  much  smaller.  In  addition,  we  observe  that 
although  the  energy  release  rate  shows  a  monotonic  dependence  on  X,  the  dependence  on 
p,  for  small  p,  is  not  monotonic. 


4.  Conclusions 

An  analytical  form  for  easy  estimation  of  the  dynamic  crack  energy  release  rate,  in 
terms  of  the  crack  tip  stress  intensity  factors,  has  been  presented  for  cracks  propagating 
along  any  of  the  principal  axis  in  orthotropic  material  systems.  The  dynamic  crack  energy 
release  rate  depends  on  the  magnitude  of  the  two  orthotropic  parameters  and  the 
instantaneous  crack  tip  velocity.  The  dynamic  crack  energy  release  rate 
is  0[(v^  -  v)  ’]  as  V  — >  ,  and  where  Vr  is  the  Rayleigh  wave  speed.  The  variation  in 

the  dynamic  crack  energy  relase  rate  for  orthotropic  materials  can  be  rationalized  in  terms 
of  the  two  orthotropy  parameters  X  and  p. 
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Figure  Caption 


Figure  1:  The  variation  of  L'\,i  (Fig. la)  and  L'\,2  (Fig. lb)  normalized  by  their 
corresponding  static  values  as  a  function  of  the  crack  tip  velocity,  normalized  by  the 
Rayleigh  wave  velocity,  is  shown  for  different  values  of  X,  and  p. 
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Table  I 


■ 

Cll 

GPa 

C22 

GPa 

C33 

GPa 

C44 

GPa 

C55 

GPa 
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GPa 

C12 
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GPa 
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5 
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5 

30 

5 

5 

5 
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wm 

m 
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30 

30 

5 

5 

5 

5 

1.0 

0.58 

0.864 

Carbon/Epoxy 

0°/90® 

2a 

10 

WM 

90 

10 

10 

5 

5 

30 

9.0 

1.29 

0.854 

2b 

90 

IB 

5 

30 

5 

0.11 

1.29 

0.946 

Carbon/Epoxy 

Quasi-isolropic 

3a 

10 

10 

5 

5 

20 

3b 
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60 

MM 
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5 

EB 

0.695 

C-SiC  ceramic 
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4a 

20 
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120 

8 

■af 

0.971 

4b 

120 

8 

20 

8 

0.17 

2.49 

0.982 
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Isotropic  Comp. 

5 

80 

80 

80 

30 

30 

30 

20 

20 

20 

1.00 

1.00 

0.911 
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Figure  1;  The  variation  of  L*\,i  normalized  by  the  corresponding  static  value  as  a 
function  of  the  crack  tip  velocity,  normalized  by  the  Rayleigh  wave  velocity,  is  shown 
for  different  values  of  X  and  p. 
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variation  of  L'^2,2  normalized  by  the  corresponding  static  value  as  a 
function  of  the  crack  tip  velocity,  normalized  by  the  Rayleigh  wave  velocity,  is  shown 
for  different  values  of  X  and  p. 
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ABSTRACT 

Bridged  crack  models  using  beam  theory  formulation  have  proved  to  be  effective  in  the  modeling  of 
quasistatic  delamination  crack  growth  in  through  thickness  reinforced  structures.  In  this  paper,  we 
model  dynamic  crack  propagation  in  these  structures  with  the  beam  theory  formulation.  Steady  state 
crack  propagation  characteristics  unique  to  the  dynamic  case  are  first  identified.  Dynamic  crack 
propagation  and  the  energetics  of  steady  state  dynamic  crack  growth  for  a  Double  Cantilever  beam 
(DCB)  configuration  loaded  with  a  flying  wedge  is  next  examined.  We  find  that  steady  state  crack 
growth  is  attainable  for  this  loading  configuration  provided  certain  conditions  are  satisfied. 


KEYWORDS 

Dynamic,  Delamination,  Crack,  Bridging,  DCB,  Stitching,  Energy  Release  Rate 


INTRODUCTION 

Through  thickness  reinforcement  of  various  kinds,  including  stitched  or  woven  continuous  fiber 
tows  and  metallic  or  fibrous  short  rods,  has  been  developed  to  address  the  delamination  problem  in 
structural  composite  laminates.  Substantial  experimental  evidence  shows  that  through  thickness 
reinforcement  dramatically  alters  the  delamination  characteristics  for  the  better  under  both  static  and 
dynamic  loading  conditions.  For  static  loading,  a  fundamental  theory  based  on  observations  of 
essential  mechanisms  is  now  mostly  in  place  [1-6].  The  mechanics  of  crack  bridging  by  the  through 
thickness  tows  has  been  mapped  out,  with  governing  length  scales  and  material  parameters  identified 
[1-6].  However  equivalent  fundamental  knowledge  and  models  for  dynamic  delamination  do  not 
exist. 
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This  paper  deals  with  the  delamination  mechanics  for  through  thickness  reinforced  structures  under 
dynamic  crack  propagation  conditions.  A  beam  theory  formulation  is  adopted  and  certain  crack 
propagation  characteristics  are  identified  for  mode  I  conditions.  In  the  next  section,  we  examine  the 
energetics  of  crack  growth  for  a  through  thickness  reinforced  DCB  specimen  loaded  by  a  flying 
wedge.  The  double  cantilever  beam  (DCB)  specimen  loaded  dynamically  by  a  flying  wedge  offers  a 
relatively  simple  experimental  approach  to  analyzing  the  mode  I  dynamic  delamination  problem. 
Regions  of  stable  crack  growth  as  a  function  of  the  material  properties  of  the  through  thickness 
reinforcement,  the  size  of  the  DCB  specimen  and  the  velocity  of  the  wedge  have  been  identified. 


Beam  Theory  Formulation  and  Solution  Characteristics; 

For  a  beam  element,  the  equations  of  motion  are; 


dN  dhi 

dx  dr 

(la) 

pBh 

^  dt^ 

(lb) 

^  -0= 
dx  ~  ^  dt^ 

(Ic) 

where  u(x,t)  and  w(x,t)  are  the  in-plane  and  transverse  displacements  of  the  neutral  plane 
respectively,  is  the  clockwise  rotation  of  the  cross-section,  t  is  the  time  variable,  A  is  the  axial 
force,  0  is  the  shear  force.  Mis  the  bending  moment,  2h  is  the  total  thickness  of  the  DCB  specimen, 
B  is  the  width  of  the  specimen,  p  is  the  density,  /  (=  Bh^/12)  is  the  moment  of  inertia  and  p(w,t)  is 
the  bridging  traction  corresponding  to  the  opening  mode.  In  this  work,  the  time  dependent  bridging 
traction  p  corresponding  to  the  opening  mode  is  assumed  to  depend  only  on  the  transverse 
displacement  w.  In  the  absence  of  an  axial  force  N,  u  =  0. 


For  a  Timoshenko  beam,  the  equations  for  steady  state  motion  can  be  reduced  to  [7]: 


dX^ 


12 R  d^w 


dX‘ 


12  R 


(R-ci)(l-ci)  h 

_  1  1  d^p(w,X)  _ 

(R-cf)Eh  dX^  ^  (R-cf)(l-cf)Eh 

M  _  (P-cf)  d^w 

dX  REh  R  QX^ 


p(w,X)  =  0 


(2a) 


(2b) 


where  X  =  x-vt ,  c]  ■=  pv^ !  E ,  =  x:  G  /£ ,  v  is  the  (constant)  steady  state  velocity,  G  and  E  are 

the  shear  modulus  and  the  Young’s  modulus  of  the  laminate  and  the  dimensionless  shear  coefficient 

x=5/6  for  a  beam  with  rectangular  cross-section.  For  steady  state  dynamic  delamination,  the  velocity 
V  is  the  delamination  crack  tip  velocity. 

For  an  Euler-Bernoulli  (E-B)  beam,  where  both  shear  deformation  and  rotational  inertia  are  ignored, 
the  equation  for  steady  state  motion  reduces  to  a  simple  form  given  by; 
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(3) 


dX^ 


12  cf  d^w  12 
dX^  Eh 


jp(w,X)  =  0 


Let  us  now  consider  a  linear  bridging  law  of  the  following  type  to  represent  the  bridging  action  of 
the  through  thickness  reinforcement: 

p(w,X)  ^  Po  + (4) 
The  linear  law  particularizes  to  the  Dugdale  law  p  =  Po  (for  and  to  the  proportional  linear 

law  p  =  (for  Pq  =0).  In  the  results  that  follow,  we  non-dimensionalize  the  variables  by  the 
laminate  thickness  h(w  =  hW ,  u  =  hU  and  X  =h^).  Thus,  the  transverse  displacement  obeys: 


d-'W  ^7  d^W 


+  b^W  +  d-  =0 


(5) 


For  Timoshenko  beam: 


12  cf  R 


(l-cf)(R-cf)  (R-cf)E 


b  = 


R 


npjh 


Ul-cf)(R-cf)  E 


d  = 


R 


12  Po 


\(l-cf)(R-cf)  E 
For  Euler-Benioulfi  (E-B)  beam: 

P  =  7[i2^  ; 

The  general  solution  to  Eqn.  5  is: 


b  = 


\\ip,h 


(6a) 


(6b) 


d^  -  J-  —  -  4 


-PI- 


+  K,e 


-J-pi.  .CZTi^ 

^  2  2  +K,e  «  ^  2 


(7) 


-4^  4 


There  are  three  regimes  to  the  solution  behaviour  which  are  independent  of  the  boundary  conditions, 
and  these  have  been  identified  below  ( Note:  S  fish/ (12  kG)  ): 

•  Case  1 :  <  0  and  P'*  >  4b^  =>  Exponential  behavior 

For  Timoshenko  beam,  this  is  true  provided: 

•2  c 

and  3  (5(1  -cf)-cff>  5(1  -  cf  )(/?  -  cf ) 


py 


E  1  +  5 

For  the  E-B  beam,  the  above  condition  is  never  satisfied. 

Case  2:  P^  >  0  and  P'*  >  4b^  =>  Oscillatory  and  non-decaying  behavior 
For  Timoshenko  beam,  this  is  true  provided: 

p\>^ /E>S/(1  +  S)  and  3(5(l-c/)-c/)^  >5(l-c/)(i?-c,^) 

For  the  E-B  beam,  this  condition  is  satisfied  when:  pv^  /E>  2^S  R 

Case  3:  p'*  <  4b^  =>  Oscillatory  with  exponential  decay  behavior 

For  the  Timoshenko  beam,  this  is  true  when: 
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(8a) 


(8b. 1) 
(8b.2) 


(8c.  1) 
(8c.2) 


3{S(\-cf)-cfy  <S(l-cf)(R-cf) 

-  For  the  E-B  beam,  this  is  true  when:  pv^/E<2^JYR 

The  conditions  determined  above  give  us  insight  into  when  dynamic  effects  can  significantly  alter 
the  mechanisms  of  deformation  and  the  resultant  bridging  phenomena.  For  instance,  if  the  crack  tip 
velocity  exceeds  the  condition  prescribed  in  Eqn.  8a,  oscillatory  displacement  fields  will  be 
introduced  in  the  wake  of  the  crack  and  these  multiple  oscillations  could  lead  to  crack  face 
interpenetration.  When  such  oscillations  are  present,  the  mechanics  of  bridging  and  the  efficacy  of 
through  thickness  bridging  ligaments  on  the  energetics  of  crack  growth  will  be  considerably  altered. 
For  example,  stick-slip  propagation  modes  would  appear  to  be  possible,  as  contacting  fracture 
surfaces  bounce.  The  complex  details  of  such  a  possibility  will  be  considered  elsewhere.  Here  we 
model  the  arms  of  the  DCB  specimen  as  a  EB  beam  and  study  propagation  characteristics  up  to  the 
point  of  fracture  surface  contact,  which  is  a  simpler  problem.  (Constants  p,  b  and  d  are  eiven  in 

Pnn  ^ 


Figure  1:  Schematic  of  through-thickness  reinforced  DCB  specimen  loaded  with  a  flying  wedge 

The  double  cantilever  beam  (DCB)  specimen  loaded  dynamically  by  a  flying  wedge,  of  constant 
ve  ocity  V,  offers  a  relatively  simple  experimental  approach  to  studying  the  mode  I  dynamic 
e  animation  problem  (Fig.  1).  The  test  is  especially  attractive  for  studying  the  bridging  effects 
supplied  by  through-thickness  reinforcement  (e.g.,  stitches  or  rods)  in  laminates.  In  figure  1,  2a  is 
e  we  ge  ang  e,  /  is  the  distance  between  the  wedge  and  the  crack  tip  and  ao  is  the  length  of  the 
bridging  zone.  In  non-dimensional  form,  I  =h  L,  and  sh  .  The  role  of  the  bridging  on  the 

crack  energy  release  rate  is  determined  in  this  section.  We  assume  that  the  crack  propagates  under 
steady  state  conditions  and  confirm  the  possibility  of  steady  state  propagation  by  finding  consisent 
so  utions.  urther,  we  assume  the  bridging  zone  size  is  invariant  and  translates  with  the  crack  tip. 
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For  the  unbridged  portion,  the  deflection  profile  =hW^)  is  obtained  by  setting  b  =  d  =  0. 
Therefore: 


^  ^  =0  for(-Z^^<-^,)  (9a) 

;^-lV 

^  ^  +  b^W  +  d^=0  for(-A,<4^0)  (9b) 

og  og 

The  relevant  boundary  conditions  are: 

w(g  =  o)=o,fV(g  =  o)  =  o,  w^cg  =  -l)  =  -a,  fv;(g  =  -l)  =  o.  (lo) 

The  governing  equation  (9)  together  with  the  boundary  conditions  (10)  and  the  continuity  conditions 
at  the  end  of  the  bridging  zone  (g  =  -Aq)  will  determine  the  deflection  profile  of  the  beam.  Note 
that  the  bridging  zone  length  (Ao)  will  be  dictated  by  the  critical  crack  opening  displacement 
s )  required  for  failure  of  the  bridging  ligament.  The  crack  energy  release  rate  ( ),  as 
determined  through  the  total  energy  balance  is: 


:  _  l{dUe.,  dU,  dU,^ 

Total  da  da  da  J 


(11) 


where  Uext  is  the  work  done  by  the  applied  load.  Us  is  the  strain  energy,  Uk  is  the  kinetic  energy,  B  is 
the  uniform  width  of  the  DCB  specimen,  and  a  is  the  crack  length.  For  steady  state  crack  extension 
a  =  vt,  where  v  is  the  crack  velocity  and  t  is  time.  For  the  DCB  specimen  loaded  with  a  flying 
wedge  this  reduces  to: 


^Total 


a  Eh  d^W^ 

6  dg^ 


phv^a^ 


(12) 


In  addition,  by  application  of  the  dynamic  J-integral,  the  energy  released  at  the  crack  tip  is  related  to 
the  bending  moment  Mby  [8]: 


^Tip  - 


{.oJ  j2 


d-'W 


dg^ 


4-oJ 


(13) 


In  general,  due  to  crack  bridging,  zlG  =  -  Gj-jp  ^  0.  The  energy  difference,  AG ,  is  the  extra 

work  done  in  fracturing  the  bridging  ligaments.  For  the  linear  bridging  law,  we  find  that  the 
shielding  contribution  AG  fi'om  the  bridging  ligaments  is: 

0 


AG  -  Gj^fgi  Gfjp  -  2  j  p  dxi 


(14) 


-w. 


This  result  is  identical  to  what  we  obtain  in  the  quasi-static  case  for  small  scale  bridging  conditions. 
Since  there  is  no  rate  dependence  to  the  bridging  law,  it  is  not  surprising  that  the  small  scale  bridging 
limit  relationship  is  obeyed. 


Since  we  limit  our  analysis  to  small  scale  bridging,  tow  failure  must  occur  in  the  wake  of  the  crack. 
Small  scale  bridging  is  ensured  provided  the  displacement  profile  monotonically  increases  within  trnr 
bridging  zone  from  the  crack  tip  and  the  pull-out  required  for  tow  failure  is  less  than  the  ma.x;~_~ 
crack  opening  displacement  within  the  bridging  zone.  This  condition  determines  a  criterion  for  the 
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maximum  allowable  bridging  zone  length,  A„,^  ,  which  is  obtained  by  solving 

dWr I  =  =  0 .  Therefore,  if  A,  <A„,^^,  then  fF,  ^  Kr,ncaif^  =  -A„,^)) ,  and  hence 

small  scale  bridging  condition  is  ensured. 

Detailed  calculations  of  the  deflection  profile,  the  crack  energy  release  rate  and  the  maximum 
allowable  bridging  length  can  be  computed  with  the  formalism  presented  above  for  both  the  Dugdale 
bridging  law  and  the  proportional  bridging  law.  For  instance,  when  the  bridging  ligaments  obey  the 
Dugdale  bridging  law,  steady  state  crack  propagation  with  small  scale  bridging  is  provided 
A  where  A^^  is  given  by; 

-2Z(Cos(A)-  Cos{l(l-A„^))  )  +  2  ADA„„,(Cos(A(J-A„,^))-  I) 

+  D(  Sw( X(l- 2A„,^JJ  +  Sw( XA^^)  )  +  D(  Sm( X  A„,^ )  -  Sin( X )  )  =  0 

and  where  A,=pL,  Dsd/X,  and  A„,^=A„,^/ L.  Regions  of  steady  state  stable  crack  growth 

under  small  scale  bridging  condition  can  thus  be  deduced  as  a  function  of  the  material  properties  of 
the  through  thickness  reinforcement,  the  size  of  the  DCB  specimen  and  the  velocity  of  the  wedge. 

CONCLUSIONS 

The  dynamic  delamination  cracking  behavior  and  the  energetics  of  crack  growth  in  through 
thickness  double  cantilever  beam  (DCB)  specimens  has  been  analyzed.  The  role  of  bridging  by 
stitches  or  rods  in  dynamic  crack  groAvth  was  computed  by  solving  the  bridged  crack  problem  within 
the  framework  of  beam  theory.  For  steady  state  crack  growth  conditions,  different  regimes  of  the 
solution  behavior  have  been  identified  which  would  correspond  to  different  crack  propagation 
characteristics.  Regions  of  steady  state  crack  growth  under  small  scale  bridging  condition  can  be 
deduced  as  a  fonction  of  the  material  properties  of  the  through  thickness  reinforcement,  the  size  of 
the  DCB  specimen  and  the  velocity  of  the  wedge.  This  provides  guidelines  for  design  of  experiments 

to  probe  the  efficacy  of  bridging  on  improving  the  dynamic  fracture  toughness  of  through  thickness 
reinforced  structures. 
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ABSTRACT 

A  nonlinear  fracture  mechanics  model  is  formulated  for  analysis  of  mode  I  delamination  of 
orthotropic  double  cantilever  beam  specimens  in  the  presence  of  large  scale  bridging  conditions. 
The  problem  is  solved  using  a  nonlinear  integral  equation  approach  in  terms  of  stress  intensity 
factors  at  the  crack  tip.  An  approximate  weight  function  is  proposed  and  validated  numerically 
for  a  pair  of  concentrated  forces  acting  on  the  surfaces  of  the  delamination.  The  model  accounts 
for  the  presence  of  regions  of  contact  along  the  wake  of  the  crack,  which  may  form  due  to  the 
action  of  the  bridging  mechanisms.  The  influence  of  the  orthotropy  of  the  material  on  the 
fracture  behavior  is  investigated  and  the  validity  of  approximated  solutions  based  on  beam 
theory  is  checked. 


KEYWORDS:  Nonlinear  fracture  mechanics,  delamination,  strengthening  mechanisms, 
anisotropic  material,  large  scale  bridging. 


INTRODUCTION 

Mode  I  and  inixed  mode  delamination  in  large  scale  bridging  conditions,  such  as  those  created 
by  through  thickness  reinforcement  in  composite  laminates,  shows  unusual  phenomena  of  crack 
face  closure,  crack  arrest  and  crack  propagation  with  crack  face  contact,  which  have  no 
precedent  in  the  delamination  of  conventional  tape  laminates  [1,2].  In  [3]  the  authors  considered 
a  typical  mixed  mode  geometry,  the  Mixed  Mode  Bending  specimen  proposed  by  Crews  and 
Reeder,  and  explained  these  phenomena  by  means  of  a  simple  analytical  model  based  on 
Timoshenko  beam  theory.  The  model  treats  the  delaminated  arms  of  the  specimen  as  beams  on 
an  e  astic,  generally  nonlinear,  foundation  of  Winkler  type  with  the  constitutive  laws  of  the 
springs  given  by  the  bridging  law,  which  characterizes  the  bridging  mechanism.  The  crack 
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closure  phenomenon  is  a  manifestation  of  the  oscillations  of  the  function  representing  the 
deflection  of  the  beams  in  the  wake  of  the  crack.  The  wavelength  of  the  function,  X,  sets  the 
characteristic  length  scale  of  the  problem,  which  in  the  case  of  linear  bridging  mechanisms  is 

given  approximately  by  7J4  =  n/2  ,  with  P3  the  modulus  of  the  foundation  and  the 

flexural  stiffness  of  the  beam  cross  section.  Once  the  limit  configuration  for  crack  tip  closure  is 
approached,  the  fracture  response  of  the  specimen  will  depend  on  the  geometry  and  the  loading 
conditions.  In  the  case  of  a  specimen  symmetric  about  its  midplane  and  in  the  absence  of  mode  II 
loading,  the  crack  will  stop  and  the  specimen  will  break  by  mechanisms  other  then  delamination. 
In  the  presence  of  mode  II  loading  or  in  asymmetric  specimens,  the  crack  will  continue  to 
propagate  and  the  propagation  will  be  opposed  not  only  by  the  bridging  mechanism  but  also  by 
friction  acting  in  the  regions  of  contact. 

The  model  proposed  in  [3]  explains  qualitatively  all  the  problems  associated  with  large  scale 
bridging  delamination.  However,  the  model  makes  strong  assumptions  which  could  affect  the 
solutions  quantitatively,  namely  it  schematizes  the  specimen  as  a  one-dimensional  structure,  it 
neglects  the  influence  of  the  elastic  material  in  front  of  the  crack  (built-in  ends  assumption)  and 
it  deals  only  approximately  with  regions  of  contact  between  the  delaminated  faces  and  the  effect 
these  regions  may  have  on  crack  propagation  driven  by  mode  II  loading. 

In  this  paper  a  nonlinear  fracture  mechanics  model  is  formulated  for  analysis  of  delamination 
crack  growth  which  removes  the  above  mentioned  assumptions,  assumes  a  two-dimensional 
deformation  field  and  accounts  for  the  orthotropy  of  the  material.  The  problem  is  solved  through 
an  integral  equation  approach  in  terms  of  stress  intensity  factors  at  the  crack  tip.  Since  the  crack 
closure  phenomenon  is  controlled  by  the  mode  I  response  of  the  laminate,  focus  in  this  initial 
work  is  restricted  to  the  problem  of  a  double  cantilever  beam. 


FRACTURE  PARAMETERS  IN  ORTHOTROPIC  DOUBLE  CANTILEVER  BEAMS 
Stress  intensity  factors 

An  exact  solution  for  the  stress  intensity  factor  iiTip  due  to  a  pair  of  concentrated  forces  P  applied 
per  unit  width  onto  the  crack  faces  of  a  double  cantilever  beam  at  a  distance  d  from  the  crack  tip 
has  been  obtained  by  Foote  and  Buchwald  [4].  They  solved  the  problem  by  applying  the  Wiener- 
Hopf  technique  to  an  isotropic,  arbitrarily  loaded  infinite  strip  and  representing  the  concentrated 
loads  in  terms  of  the  Dirac  delta  function.  A  simple  formula  approximating  the  exact  solution, 
which  has  an  accuracy  of  1.1%  and  can  be  applied  to  double  cantilever  beams  with  an  uncracked 
ligament  c  >  2h,  is  also  given  in  [4]: 


P 


=Vi2 


r. 


-+0.673 

j 


(dT 

0.815  -  +0.429 

UJ 


(1) 


where  h  is  the  half  thickness  of  the  specimen.  For  large  d/h,  Eq.  (1)  approaches  the  elementary 
beam  theory  solution  of  a  double  cantilever  beam  with  built-in  ends,  Kih^'^IP  =  -\/l2  d/h  [3].  For 
d/h  >  0.3  the  exact  ATip  is  well  represented  (error  always  lower  than  4%)  by  Gross  and  Srawley’s 
boundary  collocation  solution,  approximately  given  by  the  first  bracket  term  on  the  right  hand 


92 


BEAM  THEORY  AND  WEIGHT  FUNCTION  METHODS 


side.  The  same  limit  solution  is  given  by  the  modified  beam  theory  of  Kanninen,  which  removes 
the  assumption  of  built-in  ends  to  account  for  the  elasticity  of  the  uncracked  ligament  For  very 
small  d/h  the  dimensionless /Tip  ofEq.  (1)  approaches  Irwin’s  solution  for  a  semi-infinite  crack  in 
an  infinite  sheet,  K\h^^lP  =  (21%  h/df  ^  A  lower  limit  for  the  normalized  crack  length  a/h  of  the 
double  cantilever  beam  specimen  must  be  set  for  Eq.  (1)  to  be  valid  for  aW  d/h,  0  <  d/h  <  a/h 
Irwin’s  solution  for  very  small  ddi  is  correct  only  if  a  »  d,  and  should  be  replaced  by  Tada’s 
solution  [5]  for  a  finite  crack  of  length  a  in  a  semi-infinite  sheet  when  a/h  also  becomes  very 
small.  A  conservative  lower  limit  for  a/h  can  be  set  zsa/h>  0.3,  so  that  when  d/h  =  a/h  =  0.3, 
Gross  and  Srawley’s  solution  is  already  approached. 

The  stress  intensity  factor  due  to  a  pair  of  concentrated  forces  P\  applied  on  the  crack  faces  of  an 
orthotropic  double  cantilever  beam  at  the  coordinate  Xx  =  xu  (Fig.  l.a)  can  be  deduced  from  the 
expression  of  the  strain  energy  release  rate  Gi  obtained  by  Suo  et  al.  [6]  making  use  of  the 

orthotropic  relationship  K]i  =  -yjo^ .  Plane  stress  conditions  are  assumed  along  with  a 

principally  orthotropic  material  with  E[  =  (yl2^A^^^)/^lu^ ,  X  =  E^/Ei  and 

p  =  ^ExEj  /2Gi3  -  VfisVsi  the  orthotropic  ratios  [6],  Ei  and  £3  the  Young’s  moduli  in  the  xi  and 

X3  directions,  Gn  the  shear  modulus  and  V13  andvjj  Poisson’s  ratios.  The  dimensionless  stress 
intensity  factor  is  then  given  by; 


where: 


Pi 


+Y,(/j)r‘'0 


(2) 


Yi  (/?)  =  0.677 +  0.146(/7-1)- 0.01 78(/?  - 1)^  +0.00242(p-l)^  (3) 

and  n  =  ^(\  +  p)/2.  For  an  isotropic  material  (A  =  p=l),  Yi(p)  =  0.677  and  Eq.  (2)  coincides 
with  Gross  and  Srawley’s  solution.  The  last  term  on  the  right  hand  side  of  eq.  (2)  describes  the 
influence  of  the  elasticity  of  the  uncracked  ligament  ahead  of  the  crack  tip  and  it  vanishes  for 
large  (a-xi{)/h,  when  the  solution  for  an  orthortropic  beam  with  built-in  ends  is  recovered. 
Equation  (2),  as  well  as  the  expression  of  Gi  from  which  it  was  derived,  has  1%  accuracy  for  all 
(a-jrii)/Zi  >  2  X*’'Snd  0  ^  p  <  4. 

M  exact  solution  for  the  stress  intensity  factor  Kn  when  (a-xi{yh  <  2X^"'  is  not  available  in  the 
iterature.  However,  examination  of  Eqs.  (1)  and  (2)  along  with  the  observation  that  Irwin’s 
solution  for  very  small  d/h  =  (a-  xii)/h  maintains  its  validity  also  in  orthotropic  sheets  [7] 
suggests  the  following  formula: 
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which  has  the  right  asymptotic  behaviors  for  large  and  small  {a-xi\)/h.  The  validity  of  Eq.  (4)  for 
intermediate  values  of  (a-x\^  has  been  checked  through  finite  element  calculations  for  a  range 
of  X  and  p  typical  of  composite  laminates,  0.025  ^X<\  and  0  <  p<  4,  and  the  relative  error  has 
been  found  to  be  always  lower  than  4%.  The  lower  limit  for  the  normalized  crack  length  for 
which  Eq.  (4)  is  valid  can  be  defined  by  referring  to  the  limit  for  an  isotropic  material  and 
exploiting  orthotropic  rescaling  of  lengths  [6],  which  yields  a/h  >  0.3  When  a  reason  exists 
for  studying  very  small  cracks  as  well  as  non-small  cracks  in  numerical  work,  Eq.  (4)  can  be 
combined  very  easily  with  Tada’s  result  for  the  appropriate  domains  of  ath. 


''  Pl  Pj 

y 

p 

_ /p _ _ 

0  k 

2h 

t  t  Xi 

zn 

xii 


a 


(a) 


a 


(b) 


Figure  1 :  Schematic  of  the  DCB  specimen  under  different  loading  conditions. 

Crack  opening  displacement 

The  crack  opening  displacement  1/3  at  the  coordinate  X\{  due  to  a  pair  of  opening  forces  Pj  acting 
at  xij,  Fig.  1,  is  obtained  from  the  localized  compliance  lij  =  M3(xii)  /  Pj  which  can  be  defined 
through  an  energy  balance  or  Castigliano’s  theorem  as  shown  in  [8]  for  an  isotropic  body.  The 
localized  compliance  is  given  by: 


P.  F’  J  PP-  ^  ^  ^ 

-^1 0  «  j 

where  is  the  orthotropic  constant  defined  above,  Pi  is  a  pair  of  fictitious  forces  acting  at  xu, 
and  Xii  and  Kij  are  the  stress  intensity  factors  at  the  crack  tip  due  to  Pi  and  Pj,  respectively,  given 
byEq.(4). 


ORTHOTROPIC  DOUBLE  CANTILEVER  BEAM  WITH  LARGE  SCALE  BRIDGING 

The  stress  intensity  factor  at  the  crack  tip  of  a  double  cantilever  beam  with  tractions  p  acting 
along  the  bridged  portion  of  the  crack  as  shown  in  Fig.  1  .b  is  given  by: 

K,  (6) 

flo 

where  ao  is  the  unbridged  length  of  the  crack,  is  the  stress  intensity  factor  due  to  the  externa! 
loads  P  acting  at  .vu  =  0,  and  K\p  is  the  stress  intensity  factor  due  to  opening  tractions  /?;  Ki^Pi 
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represents  the  Green’s  function  of  the  problem  and  is  given  by  Eq.  (4).  The  tractions  p  depend  on 
the  crack  opening  displacement  and  are  a  priori  unknown  in  Eq.  (6).  If  ?/3(x,i)  >  0,  then/7[M3(xH)] 
is  the  closing  traction  developed  by  the  bridging  mechanisms.  The  value  of;?^  as  a 
function  of  tij  is  defined  through  the  bridging  traction  law,  psiuj),  which  is  one  of  the  data  of  the 
model.  If  W3(xii)  -  0,  then/»[«3(xii)]  =  -/7c[(xrii)]  is  the  opening  traction  depicting  the  effect  of  the 
contact  pressure.  The  contact  pressure  and  the  size  of  the  regions  of  contact  are  unknown  a  priori 
and  can  be  determined  through  a  compatibility  condition  for  the  crack  opening  displacement. 

The  crack  opening  displacement  W3(xii)  is  obtained  by  applying  the  superposition  principle  and 


which  yields: 


«3(^li)  =  «3(^li)p+«3(^li)p=4-P-p^ijM«3(^lj)]d^y 

^0 


(7) 


»3(^ii)_  K,,{alh,XyJh)h^-^ 

h  E[h  i  P  p. 
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a^fh  h] 
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^  K,,ialh,x,,  lh)h^-^  K^{alh,xy^  lh)h'^  . 
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P: 


(8) 

Note  that  the  dimensionless  Kfs  appearing  in  Eq.  (8)  and  in  the  equations  that  follow  depend 
also  on  the  orthotropic  ratios,  JL  and  p,  as  shown  in  Eq.  (4). 

Crack  propagation  in  large  scale  bridging 

At  the  onset  of  crack  propagation  the  crack  tip  stress  intensity  factor  of  Eq.  (6)  is  equal  to  the 

mtnnsic  fracture  toughness,  Ki  =  Kic,  and  the  dimensionless  critical  load  for  crack  propagation 
takes  the  form:  t'  t' s, 
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/7[w(Xii  /h)] 

L  J 

P30 

(9) 


where /73o  is  a  normalizing  value  of  the  crack  face  tractions,  p^,  given  for  instance  by  their 
maximum  value.  The  dimensionless  number  on  the  right  hand  side  of  Eq.  (9),  p^y^  IK^^,  is  a 

measure  of  the  brittleness  of  the  structure.  Recalling  the  expression  for  E\  and  that  K\c  = 

V^ic^i  >  Eq.  (9)  can  be  modified  to  allow  direct  comparison  between  isotropic  and  othotropic 
cases: 
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^_J__ 

yfn  Kif(alh)h^^ 
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Kii(alh,Xyi/h)h^-^ 

Pi 


p[,u(xii/h)] 

P30 


(10) 


The  normalized  crack  opening  displacement  at  the  generic  coordinate  is  obtained  substituting 
P  =  Pa  into  Eq.  (8): 

_  o  Per  K,,{alh)h^-^  K,,(a/h,x,, /h)h^-^  a  .P^ 
ip  P,  K,, 

K^{alh,xyjh)h°-^  K^{alh,xyjh)h^-^  //?)] 

a^th  max[x,i /A,x,. /A]  /■  P^  h  f^Q  h 

(11) 

The  statically  indeterminate  problem  defined  by  the  nonlinear  integral  equations  (8)  and  (9)  is 
solved  for  general  bridging  laws,  psiti^),  through  a  discretization.  A  self-consistent  solution  for 
the  crack  profile  is  obtained  iteratively  through  a  numerical  procedure  following  the  approach  of 
[8,9]. 

Dugdale  type  bridging  law,/73  =  />3o 

In  the  special  case  of  bridging  mechanisms  described  by  a  Dugdale  type  bridging  law,  p-^  =  /?3o, 
Eqs.  (9)  and  (8)  simplify  and  Eq.  (9)  alone  gives  the  dimensionless  critical  load  for  crack 
propagation.  Beam  theory  predicts  the  absence  of  regions  of  contact  for  this  case  and  this 
qualitative  characteristic  is  confirmed  by  the  more  accurate  calculations  of  the  integral  equation 
approach.  Setting  aside  for  this  paper  the  interesting  question  of  the  nature  of  contact  regions 
when  they  do  occur  (e.g.,  for  linear  bridging  laws),  a  detailed  assessment  is  made  here  of  the 
limitations  of  elementary  beam  theory  for  predicting  crack  propagation  in  the  presence  of  large 
scale  bridging. 

Figures  2.a,  2.b  and  2.c  show  dimensionless  diagrams  of  the  critical  load  for  crack  propagation 
as  a  function  of  the  normalized  crack  length  in  a  double  cantilever  beam  specimen  with  ao  =  0. 
As  already  noted  the  curves  are  correct  only  for  a/h  >  0.3  .  Three  different  values  of p'i  are 

considered,  as  marked.  The  curves  named  (a),  (b)  and  (c)  in  each  diagram  describe  the  response 
of  an  isotropic  material,  an  orthotropic  material  with  A=  0.1  and  p  =  3  (e.g.  a  graphite  epoxy 
laminate)  and  an  orthotropic  material  with  X  =  0.05  and  p  =  5  (e.g.  a  boron-epoxy  laminate), 
respectively.  The  dashed  curves  depict  the  elementary  beam  theory  solution  (built-in  ends, 
negligible  shear  deformations).  The  dotted  curve  in  Fig.  2b,  obtained  using  Timoshenko  beam 
theory  for  an  isotropic  material,  highlights  the  influence  of  the  shear  deformations. 
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Figure  2;  Dimensionless  critical  load  versus  normalized  crack  length  in  orthotropic  DCB 


specimens,  (a)  No  bridging,  (b)  Bridging  tractions ;?3  =  0. 1  .  (c)  Bridging  tractions pa  = 

03-ylGi^Ei/h  . 


Figure  2.a  confirms  the  validity  of  elementary  beam  theory  in  unreinforced  specimens  when  a/h 
is  sufficiently  high.  The  anisotropy  of  the  material  affects  the  response  only  for  relatively  small 
values  of  a/h.  The  influence  of  the  anisotropy  of  the  material  on  the  structural  response 
apparently  seems  to  be  more  marked  in  members  reinforced  through  the  thickness  (Figs.  2.b  and 
2.c).  In  this  case,  two  different  regimes  of  behavior  are  delineated  by  a  transition  value  of  a/h  = 

1/(1.73/?3o^Si(.£'i  /h)  ,  corresponding  to  the  point  where  all  curves  cross  each  other.  If  a/h  is 
smaller  than  the  transition  value,  the  anisotropy  of  the  material  strongly  affects  the  response  and 
the  elementary  beam  theory  solution  does  not  describe  the  actual  behavior  even  qualitatively.  For 
a/h  larger  than  the  transition  value,  all  curves  tend  to  become  parallel  with  a  common  slope 

given  by  1/2  pso^Gi^Ei/h  and  the  deviation  between  the  correct  solution  and  the  elementary 
beam  theory  solution  becomes  independent  of  the  crack  length  and  given  by  l/2Yi(p)X’^'^ 
P30ylGi,Ei/h  .  However,  the  fractional  error  is  Yip)K^\alh),  which  is  independent  of  the 

intrinsic  fracture  toughness  of  the  laminate,  Gic,  and  the  magnitude  of  the  bridging  tractions  and 
coincides  with  the  analogous  fractional  error  of  the  case  with  no  bridging.  It  depends  only  the 
crack  length  and  the  degree  of  anisotropy.  This  error  is  due  to  the  assumption  of  neglecting  the 
influence  of  the  elastic  material  ahead  of  the  crack  tip  and  could  be  removed  by  using  a  modified 
beam  theory  (Kanninen’s,  Williams’s). 


CONCLUSIONS 

^  approximate  weight  function  has  been  proposed  and  validated  numerically  for  a  pair  of  point 
orces  acting  on  the  surfaces  of  a  delamination  crack  in  a  possibly  thin  orthotropic  body.  The 
weight  function  allows  mode  I  large  scale  bridging  problems  in  beams  and  plates  to  be 
ormulated  as  integral  equations  without  the  limitations  imposed  on  accuracy  by  beam  theory 
approximations.  In  particular,  the  crack  tip  singularity  will  be  properly  represented.  The  integral 
equations  can  be  solved  using  well-known,  computationally  efficient  and  accurate  methods.  The 
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weight  function  strongly  depends  on  the  anisotropy  ratio.  This  is  a  feature  of  the  plate  or  beam 
geometry. 

The  first  application  of  the  new  weight  function  to  the  problem  of  a  large  zone  of  uniform 
bridging  tractions  (the  Dugdale  bridging  model)  shows  the  ranges  of  crack  lengths  over  which 
beam  theories  of  different  order  succeed  and  fail.  The  presence  of  large  scale  bridging  is  found 
not  to  significantly  increase  the  sensitivity  of  the  solutions  to  the  degree  of  anisotropy  with 
respect  to  the  case  with  no  bridging. 

While  elementary  beam  theory  will  always  be  correct  for  sufficiently  large  crack  lengths,  alh, 
there  is  a  regime  of  small  crack  lengths,  alh  <  2,  where  rigorous  solutions  are  required,  e.g. 
based  on  integral  equation  methods.  Moreover,  there  is  a  regime  of  practical  interest  for 
laboratory  specimens,  2  <  alh  <  10,  where  elementary  beam  theory  yields  only  qualitatively 
correct  trends  and  solutions  based  on  integral  equation  methods  or  on  modified  beam  theory  are 
required  for  quantitative  accuracy. 
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